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ABSTRACT 

We explore the sensitivity of two-point-mass gravitational niicrolensing to the extended nature of 
the source star, as well as the related sensitivity to its limb darkening. We demonstrate that the 
sensitive region, usually considered to be limited to a source-diameter-wide band along the caustic, is 
strongly expanded near cusps, most prominently along their outer axis. In the case of multi-component 
caustics, facing cusps may form a region with a non-negligible extended-source effect spanning the gap 
between them. We demonstrate that for smaller sources the size of the sensitive region extending from 
a cusp measured in units of source radii increases, scaling as the inverse cube root of the radius. We 
study the extent of different sensitivity contours and show that for a microlensed Galactic bulge giant 
the probability of encountering at least a 1% extended-source effect is higher than the probability of 
caustic crossing by 40-60% when averaged over a typical range of lens-component separations, with 
the actual value depending on the mass ratio of the components. We derive analytical expressions for 
the extended-source effect and chromaticity for a source positioned off the caustic. These formulae 
are more generally applicable to any gravitational lens with a sufficiently small source. Using exactly 
computed amplifications we test the often used linear-fold caustic approximation and show that it 
may lead to errors on the level of a few percent even in near-ideal caustic-crossing events. Finally, we 
discuss several interesting cases of observed binary and planetary microlensing events and point out 
the importance of our results for the measurement of stellar limb darkening from microlensing light 
curves. 

Subject headings: gravitational lensing — stars: atmospheres — binaries: general — planetary systems 
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1. INTRODUCTION 

Gravitational lensing by a system o f two point masses 
was f irst studied theoretically by [Schneider fc Weiij 
()1986l ). who mentioned that the model may be of some 
relevance for lensing of a distant source by a binary star 
in our Galaxy. At the same time they stressed that 
their study was not motivated by its applicability, the 
model was investigated rather as a specific example of 
a non-symmet ric lens. The appl i cabili ty of the model 
was studied bv lMao fc Paczvhskil (jl99lD . who estimated 
that 10% of events observed in microlensing surveys of 
stars in the Galactic bulge would be strong binary-lensing 
events. In addition, they were the first to demonstrate 
the possibility of detecting microlensing by a star with a 
planet, as an extreme case of the two-point-mass lens. 

Although the numbers of observed strong binary- 
lensing events are somewhat lower than the first estimate, 
tens of bin ary events have bee n already publi shed by the 
MAC HO fAlcock et all [2001 and OGLE (iJaroszvnskil 
l200llJaroszvnski et al.ll2'004l2006l : lSkowron et &\\mm 
projects. Many further unpublished events are listed 
on the microlensing alert web pages of OGLE^ and 
MOA^, bringing the estimated number of detected bi- 
nary events to 120-160. Most of these events can be 
described adequately by the point-source -I- binary-lens 
model. However, events with a sufficiently dense sam- 
pling of a crossing of the binary- lens caustic require more 
advanced modeling. They are sensitive to the non-zero 
angular size of the source star as well as its surface- 
brightness distribution, usually described by a radial 



limb-darkening function. The limb-darkeni ng sensitiv- 
ity of caustic-crossing microlensing events (jWittl [l995h 
turns them into a unique tool for studying stellar atmo- 
sphere s. While single-pqint-m ass events share the sensi- 
tivity (jHevrovskv et al.ll2000l ). the very high alignment 
necessary for the point-caustic of the lens to transit the 
disk of the source reduces their frequency in comparison 
with two-point-mass events, in which caustic crossing is 
highly probable. 

In principle, with good photometry and time sampling 
most caustic crossings can be used to measure the limb 
darkening of the source star. So far seven binary mi- 
crolensing events led to a reported limb-darkening mea- 
surement: MACHO 97 -BLG-28 ( Albrow ct al. 1999a), 
MACHO 9 7-BLG-41 (lAlbrow et all I2OOOD . MACHO 
98-SMC-l (|Afonso et all 120001). OGLE-1999-BUL-23 
(lAlbrow et all l2q01aD EROS BLG-2GG0-5 (lAn et al l 
20021 [Fields et al.ll2003D . MOA 2002-BLG-33 f Abe et aj 
20031) ■ and OGLE -2002-BLG-069 (Cassan ct al. 200i 
Kubas et al.l l2005l ). In addition, at the time of writ- 
ing there are six reported planetary-mass events, most 
of which are strongly influenced by the limb dark- 
ening of the source. These inc lude OGLE 2003- 
BLG-235 / MOA 2 003-BLG-53 jBond et al.l [2004), 
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OGLE 200 5-BLG-071 dU dalski ct_alj(2005D , OGLE 2005- 
BLG-390 (jBeaulieu et al.. .200&) . OGLE 2005-BLG-169 
(iGould et all [2006h. t he two-planet OGLE-2006-BLG- 
109 |Gaudi et al."'2008'), and the latest MOA-2007-BLG- 
192 {Bennett ct al. 2008). A significant fraction of these 
events have been analyzed using a simplified model of the 
caustic crossing, because most techniques for calculating 
the exact light curve of an extended source with realistic 
limb darkening arc too demanding computationally. 
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On the theoretical side, the original paper by 
ISchneider &: Weilj ()1986l ) studied in detail the lensing by 
a system of two equal masses, discussing among other as- 
pects the structure of the caustic and the critical curve, 
point-source image positions and amplifications, the ge- 
ometry of extended-source images, and local approxima- 
tions of lensing near a caustic fold and cusp. The paper 
also includes light curves of a uniform extended source 
of various sizes cr o ssing the caustic at a fold and a cusp. 
lErdl &: Schneider! ([1993') studied lensing by two unequal 
masses at arbitrary distances along the line of sight (i.e., 
in a single lens plane or in two different planes), while 
IWitt &: PettersI ()1993f ) explored the geometry of lens- 
ing by two unequal masses with an additional shear. A 
range of other general works on two-point-mass lensing 
has been published since , of which we mention here the 
study by iDominid (jl999[ ) of various limiting cases of bi- 
nary lensing. 

However, there have been fewer theoretical works 
studying in more detail the extended-source effects in 
two-point-mass lensing. These have been mostly lim- 
ited to studies of caustic crossing using a local ap- 
proximation of the cau stic by a linear fold. Here we 
point out the works by iRhie &: BennettI (|1999f) on fold- 
caustic microlensing of a linearly limb-darkened source, 
and the more general e xplor ations of fold c rossings by 
iGaudi fc Petterd ()2002al ) andlDomini^ (|2004f ). 

Only several papers de al specifically with r nicrol ensing 
near cusps of the caustic. ISchneider fc Weil3 ()1992f ) stud- 
ied the lensing by a generic cusp primarily for a point 
source, but presented also a sa mple amphficatio n map 
for a uniform circular source. iZakharovl dHH, 119991) 
contributed simple analytical expressions for the ampli- 
fications of indi vidua l imag es of a point source lensed by 
a generic cusp. iRhid (|2002f ) st udied point -source lensing 
near the cusps of a binary lens. iGaudi (SFl^ cttcrs (2002b) 
explored the microlensing behavior near a generic cusp 
for a point source and estimated the microlensing signa- 
ture of a small extended source positioned on the two 
main axes of the cusp. 

The aim of the present work is to explore extended- 
source effects in two-point-mass microlensing of limb- 
darkened sources more generally, without resorting to lo- 
cal caustic approximations. Such an approach has been 
taken so far only in individual cases to illustrate partic- 
ular effects, such as colo r changes around a sample bi- 
nary caustic presented bv lHan fc ParkI (|2001l ). In §[2] we 
briefly introduce the basic setup of two-point-mass mi- 
crolensing and the main relevant quantities. We describe 
the method we used for computing the microlensing am- 
plification in §|3] together with the limb-darkening model 
we adopted. We explore the sensitivity to the extended 
nature of the source in § |4| numerically in § 14.11 and an- 
alytically in § 14.21 The extent of the sensitive area is 
used in § 14.31 to estimate the probability of encounter- 
ing an extended-source effect of a given amplitude. We 
study the sensitivity to differences in the limb darkening 
of the source star in § [H In § [6] we assess the adequacy 
of the linear-fold approximation on the example of two 
observed events, followed by more general comments on 
other events in § [T] We discuss the results and future 
prospects in § [8l and summarize our main findings in 
§ ini A brief note on image amplifications by a generic 
cusp is added in the Appendix. 



2. AMPLIFICATION OF A POINT SOURCE AND AN 
EXTENDED SOURCE 

In a two-point-mass gravitational microlensing event 
the angular positions of images in the plane of the sky x 
are related to the angular source position y through the 
lens equation 



y = X - 



|x-xa|^ 



f-lB-, 



X - Xb 

|x-xb|2 



(1) 



where and fiB are the masses of the lens components 
relative to their total mass M {^a + Mb = 1): and xa 
and xs denote their locations. Positions x, y, xa, and 
xb are measured in units of the angular Einstein radius 
of the compound lens 
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where Dqs and Dql are distances from the observer to 
the source and lens, respectively, Dls = Dqs — Dql, G 
is the gravitational constant, and c is the speed of light. 

It is often advantageous to rewrite all vectorial po- 
sitions appearing in equation ^ in terms of complex 
variables -z = xi -|- i x2 and C = Yi + *y2 (|Wittl 119901 : 
IWitt fc PettersI [l993[ ). The corresponding form of the 
lens equation is 



fJ-A 



ZA 



ZB 



(3) 



where the lens component positions za and zb are de- 
fined similarly as the image position z. 

To obtain the amplification Aq of the flux from a point 
source passing behind the lens, we first compute the Ja- 
cobian J(x) = 9y/c)x analytically using the lens equa- 
tion (IT|). For a given source position y we then find n 
distinct image positions Xj, j = 1 . . . n, by solving equa- 
tion ID) numerically. The total amplification is the sum 
of the amplifications of individual images 



1 



det J(xjJ 



(4) 



The point-source amplification Aq diverges if any of 
the images Xj appear on the critical curve of the lens, 
i.e., the curve along which det J vanishes. The condition 
for the vanishing Jacobian can be expressed as 



det J = 1 - |7p = 0, 



(5) 



where 7 is the shear of the lens. In the case of a two- 
point-mass lens the dependence of the shear on the image 
position in complex notation is given by 



■ 7 



fJ.A 



fJ-B 



{z - zaY {z - zbY 



|7|e 



ill) 



(6) 



with real phase V'- Setting I7I = 1 in equation ^ yields a 
fourth-order complex polynomial in z, the roots of which 
are points on the critical curve. By repeating the pro- 
cedure with • varying fr om to 27r we obtain the full 
critical curve (|Wittlll99flD . 

Two new images appear on the critical curve when the 
source enters the caustic of the lens. The caustic is found 
by a reverse mapping of the critical curve using lens equa- 
tion ([T|) or, equivalently, equation ([3]). In the case of a 
two-point-mass lens the caustic may consist of a single 
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Fig. 1. — Contour plots of point-source amplification Aq by different two-point-mass lenses as a function of source position y. Columns 
correspond to lens-component separations d marked below, and each have a different scale. Rows correspond to mass ratios q marked at 
left. The origin of each plot is placed at the center of mass of the lens; crosses denote the positions of the point masses with the heavier 
body placed on the right, and caustics are marked in blue. Contours are plotted for values Aq = 2*^ increasing from the outer regions, with 
dark contours for k = {0.25, 0.5, 0.75, 1} and white contours for more widely spaced values k = {2, 3, 4, 5}, as shown in the color bar. 



continuous curve, or of two or three se parate components, 
depending on the l ens configuration (jSchneider &: Weifl 
[198 6: Erdl fc Schne ider 1993). 

We illustrate the dependence of the point-source am- 
plification Aq on the position of the source by the contour 
plots in Figure [TJ We present nine different lens config- 
urations, obtained by combining three lens-component 
mass ratios q = (xa/ fJ-B = {1/9, 1/3, 1} with three lens- 
component separations d = |xb — x^l = {0.65,1.0,2.05}. 
The choice of separations gives us three-component, con- 
tinuous, and two-component caustics, respectively, for all 
three mass-ratio values. In binary microlensing terminol- 
ogy, these are known as the "close" , "intermediate" , and 
"wide" binary configurations, respectively. In Figure [U 



as well as throughout the paper, we set the origin of the 
coordinate system at the center of mass of the lens and 
place the two point masses on the yi axis, with the lighter 
mass (A) positioned to the left of the heavier mass (B). 

In any of the lens configurations, contours far from 
the origin would eventually coincide with the circular 
contours of a single-point-mass lens with the same total 
mass. The point-source amplification grows asymmetri- 
cally from this asymptotic region and is divergent along 
the caustic. The increase to the divergence is continuous 
only when reaching the caustic at a cusp, elsewhere Aq 
converges from outside the caustic to a finite value and 
discontinuously jumps to oo. The highest amplification 
outside the caustic thus occurs at its cusps, with nearby 
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Fig. 2. — Limb-darkening model used in this work. Left panel: the first two terms /i(r) and /2(»') of an orthonormal basis obtained 
by principal component analysis (PCA) of Kurucz's ATLAS stellar atmosphere models (Hcyrovskv 2008). Right panel: full range of limb- 
darkening profiles I{r) described by the first two principal components following equation All profiles are normalized to unit total flux, 
2iT I{r)rdr = 1. The bold inner curve corresponds to k = 0, the outer solid curves are determined by limiting values Kpi; = —0.1620 
and Kfi = 0.0902. The inner dashed curves correspond to k values spaced from by 0.03, with positive values characterizing flatter profiles 
and negative values more peaked profiles. 



contours ending abruptly at the caustic. Note the narrow 
higher-ampHfication regions connecting the facing cusps 
of the caustic components in the outer columns. Inside 
the caustic Aq drops continuously to a high-amplification 
"plateau" with > 3 (Witt.& Mao 1995). 

For actual astrophysical sources (i.e., stars, in our con- 
text) the point-source approximation breaks down in re- 
gions where Aq varies substantially and nonlinearly on 
angular scales comparable to the source size, in partic- 
ular close to the caustic. In such situations the ampli- 
fication is affected by the extended nature and, more 
specifically, the surface-brightness distribution B{y') of 
the source. The amplification of an extended source 
is obtained by integrating the point-source amplification 
^0 weighted by B over the source area Ss a-nd dividing 
the result by the unamplified source flux. Alternatively, 
A* can also be computed as the ratio of the total flux in 
the images to the unamplified source flux. For a source 
with its center positioned at yc we thus have 



A(yc) = 



S^^A,{y,+y')B{y')d?y' 
j^^B{y')d^y' 
_ /^^i3(y[x]-ye)d^x 
!^^B{y')d?y' 



(7) 



Here Ei stands for the total area covered by the images 
and y ' is measured from the center of the source. In the 
numerator of the second expression we map the image 
position X to the corresponding position on the source 
y[x] using equation ([1]). For a source passing sufficiently 
close to the caustic of the lens, the spatial dependence 
of B is thus encoded in the temporal dependence of A^. 
In such cases the surface-brightness profile of the source 
can be measured, in principle, from the light curve of the 
lensing event. 



3. COMPUTATIONAL METHOD 

Evaluation of the extended-source amplification A, 
given by equation ([7]) with an arbitrary surface- 
brightness distribution of the source is a challenging task. 
Direct integration using the first expression in equa- 
tion ^ is very inefficient, due to the root-finding neces- 
sary for evaluating Aq for each point of the source and 
the need of special treatment of its divergence in caustic- 
crossing even ts. The inverse ra.y-shooting method and 
its variations (jSchneider fc WeiflllQS^lRattenburv et all 
I2002D . based on the second expression in equation ([7]), 
require prohibitively many evaluations of equation ([T]) 
to achieve a relative accuracy of at least 10~^. Such 
a minimum accuracy is necessary for analyzing light 
curves with good-quality photometry, and is essential 
for extracting the surface-brightness profile of the source 
from observational d ata. Com putational methods based 
on Green's theorem (|Gould fc Gauchcrcl 1997; Domini3 
Il998l [200l are efficient for uniform sources, but have 
difficulties in the more realistic case of sources with gen- 
eral limb-darkening profiles. Hybrid methods, which use 
a combination of several techniques, present a more ef- 
ficient approach. For example, Dong ct al. (2006) use a 
combination of ray shooting and Green's theorem. We 
employ an image-plane inte gration method suggested by 
iBennett fc Rhid (|1996D and lVermaag (|2000D . which com- 
bines elements of the direct and ray-shooting methods. 

For a given source-center position we first find the po- 
sitions of its images from the lens equation ([3]). The im- 
age positions then serve as starting points for a recursive 
scan-line flood-fill algorithm that checks points (pixels) in 
the neighborhood to see whether they ray-trace back to 
the source disk through equation ([3]). To obtain the nu- 
merator of the second expression in equation ([7|) , appro- 
priate image pixel areas are then summed with a weight 
given by the surface brightness of their corresponding 
point on the source. The flood-fill procedure runs until 
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image boundaries are reached in all directions. When 
evaluating the amplification for a given source position, 
we have to keep in mind that two new partial images 
appear whenever the edge of the source enters the caus- 
tic. These may be disconnected from the source-center 
images and thus could be potentially missed. In order 
to take them into account we calculate the intersections 
of the edge of the source and a polygon closely approxi- 
mating the caustic, which is obtained as described in the 
text following equation ^ . For initializing the algorithm 
we then use a source point just interior of such an inter- 
section instead of the source-center position, and keep 
track of the visited image pixels in order not to count 
any repeatedly. This method can easily yield the ampli- 
fication of an extended source with an arbitrary surface- 
brightness profile, with very few redundant inversions of 
equation ([3]) and relatively few back-traced rays that do 
not hit the source. 

We model the stellar source disk by a circle with 
angular radius 0^ — p^O^ and surface brightness de- 
scribed by a radially symmetric limb-darkening profile 
B{y') — I{r), with the radial coordinate r — |y'|/P* 
scaled to the source radius. We perform most of our com- 
putations for two different source radii = {0.002, 0.02} 
in units of the Einstein radius, roughly corresponding to 
typical bulge main sequence and red giant microlensed 
source stars, respectively. Limb darkening is usually de- 
scribed by the linear limb-darkening law, 

/(r) =/o[l-i;(l- Vl-r2)]. (8) 

However, the linear law often gives a poor approximation 
of limb-darkening profiles of model atmospheres. Such an 
inaccurate model would hinder the study of details of the 
caustic crossing. Instead, we model the limb darkening 
I(r) using the first two terms, fi{r) and f2{f)i of an or- 
thonormal basis obtained by principal component analy- 
sis (PC A) of limb-d arkening profiles of Kur ucz's ATLAS 
model atmospheres (|Hevrovskvll2003l I2008D . 

/(O =ai/iW +"2/2(0 =ai[/i(r) + «;/2(r)]. (9) 

The first basis function /i(r) already gives a weighted- 
average limb-darkening shape of the set of model pro- 
files, while the second term with parameter k provides 
the first correction for variations in the shape. The ba- 
sis functions /i and /2 are plotted in Figure [2] together 
with the range of limb-darkening profiles I{r) obtained 
by their linear combination. Such a model is much better 
suited for a realistic description of limb darkening, even 
though it has the same number of degrees of freedom as 
the linear law. 

4. SENSITIVITY TO AN EXTENDED SOURCE 

In this section we explore the sensitivity of the mi- 
crolensing amplification to the extended nature of the 
source. Our primary interest here is to find out how far 
from the caustic can the effect be felt. The variation 
of amplification with limb darkening, which is primarily 
important during caustic crossing but of secondary im- 
portance further off the caustic, is studied in §[5l In most 
of this section we therefore set k = and use only the 
first principal function for the limb darkening^. 

^ Such a profile is similar to that of various ATLAS models, for 
example the _B-band profile of a model with effective temperature 



We first compute in §|4J]maps of the sensitivity for dif- 
ferent lens configurations and illustrate its general prop- 
erties. In § 14.21 we derive an analytical estimate of the 
effect valid for any source not lying directly on the caus- 
tic. Finally, in I4.3l we compare the probability of events 
with an extended-source effect of a given amplitude with 
the probability of caustic-crossing events. 

4.1. Region of Sensitivity 

For a source centered at yc we describe the sensitivity 
by the relative amplification excess over the point-source 
amplification, 

A (.r \ ^*(yc) ~^o(yc) 
yc) = irr~\ • 

Figure[3]shows contour plots of 5^^ for the lens geometries 
used in Figurc[T]and a source star with radius = 0.02. 
As anticipated, most of the area sensitive to the resolved 
source is concentrated near the caustic curves. However, 
the |(5cx| — 0.001 and 0.01 contours depart in some cases 
substantially from the caustic, as discussed further below 
and in § ^ 

The dot-dashed contour extending from the cusps cor- 
responds to a zero excess, with the extended-source am- 
plification equal to the amplification of a point source 
positioned at its center. In the dominant positive green 
areas the amplification exceeds j4o(yc). This is to be 
expected, since for a source positioned fully outside or 
inside the caustic a large part of the limb lies closer to 
the caustic than the source center, and the point-source 
amplification increases convexly toward the caustic. 

The limited negative red areas, in which the amplifica- 
tion is lower than j4o(yc), can be seen as compact lobes of 
varying extent along the outer axes of all cusps. Further 
away from any cusp along its axis 5cx changes back to 
positive. The more subtle reason for this effect is caused 
by the concave drop of amplification in both directions 
perpendicular to the axis, which is initially stronger than 
the outward convex decrease along the axis. The off-axis 
parts of the extended-source limb thus decrease the total 
amplification below the value corresponding to its center. 
Further from the caustic the relation reverses and the 
perpendicular drop becomes weaker than the radial drop. 
On a surface plot corresponding to the point-source am- 
plification maps in Figure [T] these regions would appear 
as sharp ridges extending outward from the cusps and 
eventually merging with the surroundings. Neglecting 
the change of amplification perpendicular to the cusp 
axis can lead to misleading results, s uch as the posi- 
tive ext ended-source effect obtained by iGaudi fc PettersI 
()2002br ) along the full extent of the outer axis of a generic 
cusp. For details see § 14.21 

The smallest negative lobes can be found along the lens 
axis: for non-equal-mass binaries (g < 1) the smallest lies 
around the heavier lens-component; for an equal-mass bi- 
nary (g = 1) the two lobes closest to either component 
are the smallest. Interestingly, it is in this region that the 
increase of point-source amplification toward the caustic 
is the strongest. However, as indicated by the more cir- 
cular contours in Figure [H in this case the extending 

T^j f = 12500/4" and surface gravity log g = 2.5 or the ij-band 
profile of a model with T^f j = 8000 and log g = 1.5 (both with 
solar mctallicity and microturbulent velocity vt = 2 km s~^). 
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Fig. 3. — Extended-source effect 5ox for a source with radius p, = 0.02 and ft = limb darkening as a function of source-center position 
yc. Panels correspond to the same lens configurations as in Figure [T] Solid contours are plotted for 5cx = ±0.001, ±0.01, and ±0.1; the 
dot-dashed contour corresponds to (5ox = 0. As marked in the color bar, positive values are mapped in shades of green with black contours, 
negative in shades of red with white contours. Areas with |(5ex| < 10~ are left white. The open circle in the upper left corner of each 
panel illustrates the size of the source (note the different scales in each column). 



ridge is more blunt, nearly radially symmetric in some 
cases. The perpendicular drop soon becomes compara- 
ble with the radial drop, which then dominates and leads 
to a positive excess. 

The most interesting and unexpected feature can be 
seen in the outer columns with the compound caustics. 
The negative sensitive area connects facing cusps of the 
caustic components, extending many source radii along 
the axis from the caustic. These regions correspond to 
the narrow higher- features of Figure [T] noted earlier 
in § m The amplification excess can reach well over 1% 
in these regions and is thus observationally significant. 
For a given lens separation d the significance or strength 
is dependent on the mass ratio q. With decreasing q the 



caustic components become more separated, the regions 
become more diffuse perpendicularly, and the maximum 
|(5ox| decreases. Similarly, for a fixed q the importance 
of this region depends on d. For instance, for q — X/Z 
(middle row of Figure [S]) the connecting area splits into 
opposite lobes for d < 0.41 and for d > 3.55. Note that in 
the small-separation limit the two-point-mass lens con- 
verges to a single-point-mass lens with the same total 
mass, while the large-separation limit leads to two inde- 
pendent point-mass lenses. Many events in both regimes 
are hard to distinguish from single-point-mass microlens- 
ing. 

In the top panel of Figure |4] we present an enlarge- 
ment from Figure [31 namely, the bottom right lens with 
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Fig. 4. — Extended-source effect (JoxCyc) for a lens with d = 2.05 and q = 1/9, and sources with p« = 0.02 (top panel) and p» = 0.002 
(bottom panel). Symbols and lines have the same meaning as in Figure [3] with four additional dotted contours at ±10"^'^ and ±10~^'^. 
Inset in top panel: cut along the source track marked with the orange arrow. The time axis is parameterized with t expressed in units of 
source-radius crossing time and t = corresponding to the crossing of the caustic axis y2 = 0. Values corresponding to 5ex = 0, ±10~^, 
and ±10"^'^ contours are marked in the inset by gray lines. 



d = 2.05 and q = 1/9. Details of the structure close 
to the caustic are better visible in this blowup. On 
the inner side of the caustic we may begin to discern a 
very narrow negative region, where the diverging point- 
source amplification dominates over the non-divergent 
extended-source amplification, bringing S^x in the limit 
to —1 along the inner side of the caustic. Just outside 
the caustic, at a sufficient separation from cusps, the 
higher positive contours converge to a curve parallel to 
the caustic at a distance of one source radius. Here the 
extended-source limb enters the caustic while the point 
source remains outside, thus only the extended-source 
amplification starts to grow rapidly. The area between 
this curve and a similar parallel curve inside the caustic 
is the key region where extended-source effects have been 
considered seriously so far. From our results it is clear 
that at the 1% level this region should be substantially 
expanded close to cusps, especially along the axis con- 
necting the caustic components. In this example the re- 
gion around the smaller caustic component on the right 
is fully dominated by cusp proximity, so that the one- 
source-radius parallel curve is relevant only at the 10% 
level. At the level of 0.1% the sensitive region broadly 
covers the entire caustic structure including the area be- 
tween the components. 
In order to demonstrate the large extent of the sen- 



sitive zone we included an inset plot illustrating the Sex 
variation with time for a source traveling along the arrow 
marked in the main plot. The maximum effect reaches 
nearly —1% when the source is located along the axis 
17.6 /9* from the nearest caustic point (—1% would be 
reached 16.2 from the cusp). Hence, even in such en- 
tirely non-caustic-crossing events, neglecting extended- 
source effects may lead to inaccurate light-curve fits and 
biased event parameters. 

In the bottom panel of Figure[3]we illustrate the scaling 
of Sex with source size by using the same lens parameters 
and a ten times smaller source, p^, = 0.002. Here even 
the outer 0.1% contour closely follows the caustic. As ex- 
pected, in comparison with the 1% and higher contours 
from the top panel, the distance of the contours from the 
caustic generally scales with the source size. Neverthe- 
less, the scaling is not uniform, as can be noticed from 
the somewhat more elongate non-zero contours along the 
cusp axes. Indeed, for a source crossing the caustic axis 
a —1% effect occurs as far as 28 p* to the right of the 
cusp of the larger caustic component. 

This can be seen also in Figure O which includes a 
detail of the top center panel of Figure [3] (with q = 1 
and d = 1) and a further blow-up of Figure 21 each for 
source radii 0.02 and 0.002. In the top left panel we in- 
clude cuts through the plot along three marked source 
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Fig. 5. — Details of extended-source effect (Scx{yc) for a d = 1, 
q = 1 lens (top panels), and a d = 2.05, q = 1/9 lens (middle 
and bottom panels) for source radii p, =0.02 (top left and middle 
panels) and p» = 0.002 (top right and bottom panels). Meaning 
of symbols and lines as in Figure |4] Top left inset: cuts for source 
tracks marked by arrows; top right inset: cuts for the same number 
of source radii along the same tracks. Caustic crossing of the source 
center occurs at time t = 0; t» is the source-radius crossing time. 
Gray horizontal lines in the insets mark contour values (Sex = ±0.1. 



trajectories crossing the caustic. In the top right panel 
we inchide cuts for the smaller source along the same 
straight lines extending the same number of source radii 
from the caustic (i.e., in this panel the arrows would be 
shrunk ten times towards the caustic). The A trajecto- 
ries correspond to the well-studied fold-crossing regime. 
As the source approaches the caustic from inside, the 
positive effect peaks just after the limb touches the caus- 
tic and drops to ^ox = — 1 at the caustic. The peak 
value is only slightly higher for the smaller source. Once 
the source center is outside the caustic, the correspond- 
ing point-source amplification jumps to a low value and 
the extended-source amplification dominates by a fac- 

— 1/2 

tor roughly proportional . Nevertheless, S^y: drops 
rapidly practically to zero when the source limb fully 
exits the caustic. The region with |(5ox| ?J 0.01 is thus 
limited to a band from 2.4 inside to 1 outside the 
caustic. For the region with \Scx\ ^ 0.001 the outer limit 
remains the same, but the character of the inner region 
in this particular situation depends on the source size. 



While for the smaller source the inner limit of the band 
hes at 7.8 p*, for the larger source the width of the entire 
caustic is less than double (only 15 p*) along the hori- 
zontal line extended from A. The S^x value thus remains 
larger than 0.001 along the line inside the caustic. The 
validity of the fold-caustic approximation is discussed 
further in § [S] 

Sources moving along trajectories B directly exit the 
caustic through the cusp along its axis. The character 
inside the caustic is similar to A with a few important dif- 
ferences. The peaks are substantially higher than those 
of the corresponding A trajectories. The positive peak 
occurs already 2.4 p* from the cusp for the larger source, 
when the limb of the source touches the inclined parts of 
the caustic above and below the axis. Due to the curva- 
ture of the caustic, for the smaller source the peak occurs 
as far as 4.9 p* from the cusp. Outside the caustic the 
curves are entirely different from A, slowly rising from 
— 1 towards zero, the larger source passing 6cx — —0.1 at 
2p*, (5ex = -0.01 at 4.4 p*, and (5ox = -0.001 at 7.6 p* 
from the cusp. The rise is even slower for the smaller 
source, which passes i5cx = —0.1 at 4.6 p*, 6ex = —0.01 
at 11 p*, and i5ex = —0.001 as far as 23 p* from the cusp. 
Inside the caustic the larger source crosses (5ex = 0.01 at 
5.9 p* and i5cx — 0.001 at 14 p* from the cusp, the smaller 
source at 12 p* and 26 p*, respectively. The important re- 
sult that the smaller source size leads to a larger extent 
of the sensitive region in terms of source radii on both 
si des of the cusp contradicts one of the analytical findings 
of lGaudi fc Petteri ()2002bl ). who concluded that the ex- 
tent was independent of source size for sufficiently small 
sources. We discuss this point further in § 14.21 

Source-center trajectories C pass through the cusp per- 
pendicular to its axis. Except for the immediate vicinity 
of the cusp, the values are positive, peaking at Sex — 1 at 
positions 0.65 p* from the cusp, dropping to 0.1 at 1.2 p* 
from the cusp, to 0.01 at 3.1 p*, and to 0.001 at 8.2 p* 
from the cusp. The decay rate is thus intermediate be- 
tween cases A and B. In this case the curves for the two 
source sizes are nearly identical, differing only slightly in 
the peak height. 

The source-size dependence of the extended-source ef- 
fect is illustrated in more detail in Figure [6] for the three 
cuts from Figure [5] and a source-center-crossing single- 
point-mass lens for comparison. The graphs for cuts A, 
C, and for the single lens are plotted as a function of 
time in units of the source-radius crossing time t^. For 
cut B time is shown in units oftB(p*) = (p*/0.02)~^/'^ t*, 
in order to account for the changing characteristic scale 
along the cusp axis discussed above. In each panel the 
curves with increasing peak height correspond to de- 
creasing source size. For cut A, the values of p* = 
{0.1,0.05,0.02,0.01,0.005,0.002,0.001}, with the previ- 
ous caustic entrance visible in the inset for p* — 0.1; for 
cuts B and C, p* = {0.1,0.05,0.02,0.01,0.005,0.002}; 
for the single lens, p* = {2,1,0.5,0.2}. The Sex curves 
converge for p* ^ nearly everywhere to limit curves, 
so that curves for smaller source radii would be indistin- 
guishable in the plots. The only obvious exception oc- 
curs for a source centered within one radius outside the 
caustic in cut A, where Ao(t/t*) converges to a constant 

— 1/2 

while A^{t/t^) and thus also Scx{t/t*) diverge as p* . 
It is interesting to note the rapidity of the convergence in 
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Fig. 6. — Source-size dependence of extended-source effect <5ex- Panels correspond to the source tracks marked in Figure|5]as noted in 
the top left corner, and to a source-center-crossing single-point-mass lens (bottom right) for comparison. Panel A inset shows an expanded 
view of the plot. Curves with increasing peak height correspond to source radii p, = {0.1,0.05,0.02,0.01,0.005,0.002,0.001} for panel 
A, p, = {0.1,0.05,0.02,0.01,0.005,0.002} for panels B and C, and p. = {2,1,0.5,0.2} for the single lens. Time is marked in units of 
p*-crossing time tt except in panel B, where it is marked in units of the rescaled teCp*) = (p*/0.02)~^''^ 4* to account for the changing 
characteristic scale along the cusp axis. Vertical dotted lines mark the times of first contact of the source with the crossing point of the 
caustic, source-center crossing, and last contact. In panel B the first and last contact with the cusp occurs at t/t^{p*) = =f(p*/0-02)^/''. 
Horizontal dotted lines mark <5cx = ±0.1. 



the case of the single-point-mass lens, in which the hmit 
curve i s an a nalog of the uniform-source B{z) curve of 
iGouldl ([1991 . In the case of the three presented binary- 
lens cuts the limit curves are valid for sources at least 
two orders of magnitude smaller than in the single-lens 
case. For instance, along the fold-crossing cut A, the 
t < Q part of the (5ex curve is independent of the source 
radius only for < 0.001, and the shape of the t < 
part of the corresponding light curve is independent 
of only for < 0.005 (not shown here^). 

Even though we explored the source-size dependence 
only for three specific caustic crossings of a specific lens, 
due to the universality of the local lensing behavior the 
cut A limit curve is valid for any fold crossing, provided 
we replace by the source-radius cawstic-crossing time 
in non-perpendicular crossings. Similarly, the cut B and 
C limit curves are valid for any parallel and perpendic- 
ular cusp crossings, respectively. However, in all cases 
the rate of convergence as — > does depend on the 
geometry of the particular crossing, primarily on the lo- 

* The faster convergence of the light-curve shape is caused by 
the partial compensation of the steeper drop of (5ex at t — > 0~ for 
a smaller source by the steeper rise of Aq closer to the caustic. 



cal curvature of the caustic, cusp shape, and vicinity of 
(other) cusps. 

The opposite large-source regime is illustrated in the 
middle panel of Figure [51 where the source size is compa- 
rable with the size of the caustic component. Hence, the 
proximity of all cusps has an influence on any caustic- 
crossing event, as demonstrated by the smooth deformed 
contours, very different from those in the top left panel. 
Along the lens axis the 6cx — —0.01 contour extends 
14.6/9* out from the leftmost cusp, but only 1.29 p* out 
from the rightmost cusp. The structure for the smaller 
source in the bottom panel looks more familiar. Here the 
ratio between the (5ox = —0.01 left and right contour ex- 
tent is relatively less extreme but still large: 26.2 from 
the leftmost cusp and 4.8 p* from the rightmost cusp. 
Note that the only significant region inside the caustic 
with a negligible extended-source effect \6ex\ < 0.001 lies 
within the innermost contour. 

To summarize the influence of source size, a smaller 
source reaches higher ampliflcation close to the caustic. 
The amplification excess |i5ex| typically also increases, but 
for sufficiently small sources on the inner side of folds and 
in any direction from cusps it remains constant. While 
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the extended-source effect for a smaller source is rele- 
vant in smaller regions around the caustic in terms of 
Einstein radii, it extends substantially further outwards 
from cusps in terms of source radii p,. Along the axis 
connecting caustic components a 1% effect may occur 
tens of from cusps for the smaller source, and signifi- 
cantly over 10 even for the larger source. 

4.2. Analytical Estimate of Extended- Source Effect 
4.2.1. General Results 

For any sufficiently small source we may study the am- 
plification excess 5cx analytically by expanding the point- 
source amplification appearing in equation ([7]) around . 
In principle, such an approach should give us good results 
for any region where Aq (y) is analytic within a source ra- 
dius p, of point y^ , which means anywhere except within 
a source radius of the caustic. The extended-source am- 
plification can then be written as 



0.5 



i=i 

1 ^ d^A^ Ly'^y',B{y')d^y' 



+ 0{pl). ill) 



The first correction to the point-source amplification is 
given by the dot product of the point-source amplifica- 
tion gradient and the position of the center of brightness 
of the source measured from source center, 



J^^y'Biy')d'y' 
/j.^S(y')d2y' 



Yb 



(12) 



where yb is the position of the center of brightness mea- 
sured from the origin of the global coordinates. The fol- 
lowing terms in equation pip involve higher-order deriva- 
tives of the point-source amplification and higher-order 
moments of the brightness distribution. 

In the special case of a circularly symmetric brightness 
distribution B(y') = I{r), which is of primary interest 
here, we have yb = yc and all terms of odd order drop 
out and the first three terms of the resulting expansion 
are 

A(yc) = ^o(yc) + - AAo(yc)p^ \ \ , + 

4 I{r)rdr 



1 

64 
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I{r)rdr 



+ 0(pS),(13) 



where A is the Laplacian and A^ = AA is the bihar- 
monic operator. For a symmetric source the leading- 
order approximation of equation pop therefore is 

4Ao(yc) /(r)rdr 

This expression provides an easy way to estimate the sen- 
sitivity to an extended source merely from the knowledge 
of the point-source amplification^. 

^ Its Laplacian can be easily computed numerically from a map 
of A)(y) such as those in Figure [T] 




Fig. 7. — Comparison of exten ded- source effect <5ox(yc) and its 
approximation given by equation H14| l for the lens and source from 
the top left panel of Figure [5] The caustic is marked by the bold 
blue line, exact contours are plotted in black and approximate con- 
tours in red; contour levels are the same as in Figure |5] 

Except for the immediate vicinity of the caustic, the 
relative amplification excess can thus be approximated 
by a product of three terms, each depending on different 
parameters. Only the first term depends on the source- 
center position and the lens parameters, therefore it fully 
determines the contour geometry. The contour shape 
away from the caustic is thus independent of the source 
size and limb darkening. The second term defines the 
quadratic scaling with source radius, i.e., doubling the 
source size at a given position will increase 5cx four times. 
The last term depends purely on the limb darkening of 
the source and varies only weakly. For the linear law 
given by equation ([5]) the term equals (1.5 — 0.7v)/{3 — 
v), for our PCA model given by equation ([9]) we get 
(0.2155 -|-0.3097k)/(0.4665-|-0.3272k). For either model 
these values range from 0.4 to 0.5 from the most peaked 
to the flattest possible profile. Such a weak dependence 
on the limb darkening away from the caustic justifies our 
usage of a fixed limb-darkening profile in this section. Its 
parameter k = gives us 0.462 as the value of the last 
term. For orientation, the most peaked profile would 
have \6cx\ lower by 13%, the fiattest would have \Scx\ 
higher by 8% away from the caustic. 

In Figure [7] we compare this approximation to the true 
(5cx using the detail in the top left panel of Figure [5] The 
agreement is remarkable anywhere more than a source 
diameter from the caustic and very good up to a source 
radius from the caustic except near the cusp. In the 
outer vicinity of the caustic the exact contours bend and 
run parallel to the caustic at a distance of p*, while the 
approximate contours extend directly to the caustic as if 
Sex continued analytically. The slight discrepancy closer 
to the cusp corresponds to the change of contour geome- 
try with source size in this region, as noted earlier in the 
discussion of Figures [4H6l Adding the higher-order term 
from equation (jl3p would improve the agreement here^. 
Finally, note that all outer approximate contours extend 

^ This wa s demonstrated after the manuscript submission by 
IGouldl II2008I ). whose hexadecapole approximation is a numerical 
implementation of equation II13I I including the 0{pf ) term. 
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directly to the cusp, while all exact contours extend to 
points slightly offset to both sides of the cusp. 

The zero-effect contour clearly has a special signifi- 
cance. Equation (|14p shows that away from the caus- 
tic Sex — corresponds to points with Ay4o(yc) = 0, 
independent of the source properties. This explains the 
similarity of the zero contours for both values of in 
Figure [4] or Figure [Sj In fact, the similarity allowed us to 
use the zero-Laplacian curve instead of the zero-<5ox curve 
in several of the previous plots further from the caustic, 
as it is less sensitive to low-level numerical noise. We also 
point out specifically that away from the caustic the zero- 
effect contour remains the same for any limb-darkening 
profile I{r), as discussed further in §[5l 

For completeness, we briefly turn to the case of an 
asymmetric brightness profile, relevant in our case for ex- 
ample for a spotted source star. We can simplify equa- 
tion (jlip similarly by expanding Aq around the center 
of brightness yb instead of the source center yc. This 
approach again leads to the elimination of the first cor- 
rection and we are left with 



7i:„ y;yLB(y')dV' 



Yhcj Ybcfc 



+ Oipi), (15) 



where ybc is given by equation (|12p . The amplification 
of an asymmetrically bright extended source lying off the 
caustic is thus equal to the amplification of a point source 
placed at its center of brightness plus a correction of 
second order in source radius. 

We emphasize here that the expansions in equa- 
tion ()lip and equation ()15p are valid for any gravita- 
tional lens and a small source of any shape and bright- 
ness distribution, not just a binary lens and a circular 
source. Similarly, equation (I13p and equation ()14p are 
valid for any lens and any small source with a circularly 
symmetric brightness distribution. For instance, in the 
special case of a single-point-mass lens and a circular or 
elliptical source equation (fTTI) yields the analytical result 
obtained by Hevrovsky & Loeb (1997). Equation (|13p 
applied to a single-point-mass lens and a circular sym- 
metric source re p roduc es the analytical result given by 
iHevrovskv eTall ()2000f ). 

4.2.2. Extended Source Near a Cusp 

In order to examine the peculiarities of extended- 
source effects near cusps, we start from the analytical 
expressions for the amplification of im ages formed by a 
gener ic cusp, in the form derived by IZakharov] (jl995l . 
Il999f ) and presented here in the Appendix. For a suf- 
ficiently small source lying off the cusp we may then 
utilize equation ([T4|) derived above to obtain i5ox- We 
express the results here in terms of the source position 
(yiliYJ-) measured in a local coordinate system centered 
on a given cusp, with the first axis identified with the 
cusp axis and the second perpendicular to it. As noted 
in the Appendix, without loss of generality we rotate the 
coordinates so that the cusp is pointed towards negative 
y|| and opened towards positive y|| . 

We are primarily interested in getting expressions for 
(5cx of a small source placed along the cusp axis {y± = 0) 



outside (y|| < 0) or inside (y|| > 0) the cusp, or perpen- 
dicular to the cusp axis (with y|| = 0). For this purpose 
we compute the Laplacian of the total amplification , 
express Sex using equation (I14p . and get our results by 
setting y^ = or yy — 0, as required. 

Outside the caustic the generic cusp approxim atio n 
yields a single image and is given by equation (IA7p . 
Starting with the most interesting region along the outer 
cusp axis, we arrive at the exact expression 



^tot 



, 8K 2 
y± = 0,y|| <0 =-3 + — , 
9 y^ y| 



(16) 



where the cusp parameter K defined in equation (jA3p is 
always positive in our coordinate system. Note that the 
first term, which dominates in the vicinity of the cusp 
and originates from the second derivative in the perpen- 
dicular direction, is negative, because yy < along the 
outer cusp axis. From the leading order term we readily 
get 



2 K 

Scx{y± = 0) « 77-3 p: 



^y\\ Ionr)rdr ' 



(17) 



a negative extended-source effect along the axis dropping 
off with the inverse cube of the distance from the cusp. 
The effect is directly proportional to the cusp parameter 
hence i5ex is suppressed near broad cusps with small 
K (such as the cusp close to the heavier component of the 
binary) and boosted near narrower cusps with large K , 
which explains the large extent of the sensitive zone be- 
tween the narrow facing cusps of multi-component caus- 
tics. Expressing the distance in source-radius units, we 
see that at a given number of source radii from the cusp 
(5cx is inversely proportional to the source radius . By 
inverting equation (|17p we get the scaled contour extent 



^(y. 



0) 



2K ^ll{r)r^dr 



-I 1/3 



9p*(5ox I{r)rAr 



(18) 



Hence, even though for smaller sources the extent of a 

given ^ex contour along the outer cusp axis shrinks as 

2/3 1/3 

pj when measured in Einstein radii, it grows as 

when measured in source radii, explaining our observa- 
tions from § 14.11 Measured either way, the extent also 

scales as Ocx 

The expressions derived above are in good agreement 
with the numerical results presented in § 14.11 Their 
disagreemen t with the r esults obtained along the outer 
cusp axis bv lGaudi fcTe ttcrs (2002b) and Gould (20081) 
comes from our accounting for the change of the ampli- 
fication in the direction perpendicular to the axis in its 
immediate vicinity. As mentioned earlier in § 14. li close 
to the cusp this change dominates over the relatively slow 
drop al ong t he axis . The approach of Gaudi & Petterj 
(2002b) and iGouldl (j2008|) takes into account only the 
radial drop and thus yields only the smaller second term 
in equation (I16p . which in turn incorrectly implies a pos- 
itive extended-source effect that is independent of cusp 
properties. The obtained contour-extent scaling would 
be the same as at a fold, i.e., linear with source radius 
if measured in units of Einstein radii, or independent of 
source radius if measured in units of source radii. 
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Staying outside the cusp and studying the effects of a 
small source displaced from the cusp perpendicularly to 
its axis, we arrive at the exact result 



(yii - 0) 
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For the extended-source effect we get 



5c^{y\\ = 0) 



5 2 /o^ /(r)r3dr 
18 yi ^* /o /(r)rdr ' 



(19) 



(20) 



a positive effect dropping off as the inverse square of the 
perp endicular separation frorn the cusp, in agreement 
with iGaudi fc PettersI ()2002hf ). Expressing the scaled 
perpendicular contour extent 



— (yii=o) 
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/o /(r)r3dr 
cx /(r)rdr 



1/2 



(21) 



we see it is independent of the source radius, in agreement 
with our C cuts in Figure [3 Unlike the inverse-cube-root 
scaling with (5ex along the axis, perpendicularly from the 

cusp the extent scales as Sc-^^"^ . 

Inside the cusp there are three images with total am- 
plification given by equation (|A5p . Repeating a similar 
procedure as above, we get exactly 
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(22) 



the same as the result in equation (|16p for the outer axis. 
However, along the inner part of the cusp axis y|| > 0, 
thus both terms a re positive. The corresponding result of 
iGaudi fc Petterd ([2002b) includes only the smaller sec- 
ond term. 

Not surprisingly, equations (fTT)) and (fT8)) are valid even 
inside the cusp, the only difference being the positive val- 
ues of 5cx- Sources positioned along the axis symmetri- 
cally on opposite sides of a cusp thus have an opposite 
value of (5ex, provided they are far enough to lie off the 
caustic. The required distance can be estimated by com- 
puting the position on the axis of a small source tangent 
to the cusp from inside. 



ytan 



K 
p* 



-,1/3 



(23) 



The tangent position of the source thus also scales as 

when expressed in source-radius units. This ex- 
plains another observation from ji l4.11 namely, the change 
with radius of the position of the first peak along cut B 
in Figure [5l Even though the peak doesn't lie exactly at 
the tangent point, it lies only slightly closer to the cusp 
where a larger part of the limb is aligned with the caus- 
tic. We note that the requirement y|| > ytan implies that 
inside the cusp equation holds for values <5cx ^ 0.08. 

Usage of the formulae derived above in actual gravita- 
tional lensing scenarios requires more than just satisfying 
the conditions of equation (fT4|) . i.e., having a small source 
lying off the caustic. Additional constraints are placed 
by the applicability of the generic isolated cusp model, 
primarily putting limitations on the distance from the 
cusp. For instance, in binary lensing there are three im- 
ages of a source positioned outside the caustic and five of 
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Fig. 8. — Projected widths lw{a) in Einstein radii of regions with 
absolute extended-source effect |<5ox| > w for source trajectories ori- 
ented at an angle a from the lens axis. Lens geometries correspond 
to the middle row of Figure = 1/3, d marked in each panel) 
with source radius p, = 0.02. Curves for w = 0.001, 0.01, and 0.1 
are marked by dotted, short-dashed, and long-dashed lines, respec- 
tively. Solid lines mark the projected width of the caustic-crossing 
region lcc{o)- The symmetry around a = 7r/2 arises from the axial 
symmetry of the lens. 

a source inside, as opposed to one and three images, re- 
spectively, in generic cusp lensing. At a greater distance 
from the cusp the additional images cannot be neglected. 
The case of facing cusps also places obvious distance lim- 
its. 

4.3. Probability of Extended-Source Effect 

Having shown that the region sensitive to the source 
size may be considerably larger than the region occupied 
by the caustic, we now quantify the increase in detec- 
tion probability. As in all similar cases, rigorous estima- 
tion of the probability of detecting an extended-source 
effect would require a light curve by light curve anal- 
ysis, based on measuring the residuals from a best-fit 
point-source binary-lens light curve. Nevertheless, un- 
derstanding that the actual probability will tend to be 
lower, we may estimate the sought probability increase 
simply by computing the number of random straight-line 
source trajectories that cross a given |(5ex| contour. For 
trajectories oriented at a given angle a from the lens axis 
(yi), the number is proportional to the perpendicularly 
projected extent of the contour. Specifically, wc define 
lw{<^) as the total projected width^ perpendicular to tra- 
jectory orientation a of the contour \Sox\ = w. Similarly, 
to estimate the frequency of caustic-crossing events, we 
define /cc(ck) as the total projected width perpendicular 
to trajectory orientation a of the outer parallel curve of 
the caustic at a distance p^. 

Sample Iwict) and Icdcn) functions are displayed in Fig- 
ure [8] for lenses with g = 1/3 and source radius p, — 0.02 
(middle row of Figure [3|) . Clearly, the greatest increase 
over caustic-crossing events occurs when contours con- 
necting multiple caustic parts are crossed perpendicu- 
larly to their maximum elongation. For instance, for hor- 

^ Sum of projected widths of all components of the contour with 
any possible regions of overlap counted only once 
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TABLE 1 

Probability Pu, of occurrence of extended-source effect of amplitude |<5cx| > in multiples of 

caustic-crossing probability 
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izontal (a = 0, tt) trajectories and the three-component 
d — 0.65 caustic, the 1% effect region is 2.3 times wider 
than the caustic-crossing region. For vertical (a = 7r/2) 
trajectories and the two-component d — 2.05 caustic, the 
1% region is 1.9 times wider than the caustic-crossing re- 
gion. 

Note that there is a sHght difference between the quan- 
tities used here and si milar ones found in the li terature. 
The width Ap used bvl Mao fc Paczyhskil ()1991[) and the 
equivalent Scc used by iNight et al.l ( 20071 ) are similar to 
Icc- However, unlike Zee these quantities purely measure 
the size of the caustic, thus including only trajectories 
for which the source center crosses the caustic instead 
of all caustic-crossing trajectories. For tiny sources the 
difference is negligible, but even for the = 0.02 source 
used in Figure [H Scc(7r/2) in the top panel is 16% (four 
source radii) lower than Zcc(7r/2). In addition, when the 
caustic shrinks to isolated points in the very wide c? ^ 1 
or very close d <^ I binary limits, Scc drops to zero while 
Icc converges to two or three source diameters, respec- 
tively^. 

The relative probability of achieving at least an \Scx\ = 
w extended-source effect in comparison with the prob- 
ability of caustic crossing is given by the ratio of the 
angle-averaged projected widths, 



p — 



(24) 



We present the values of for the nine lens geometries 
of Figure [3] in the left part of Table [T] for source radii 
p* = 0.02 and 0.002. The obtained values demonstrate 
that the probability of an extended-source effect can be 
much higher than the probability of caustic crossing. In 
the case of the larger source the probability ratio for a 
1% effect exceeds 2 in two of the nine studied geometries. 
For the smaller source Pq.oi > 1-3 also in two of the nine. 
Apart from the obvious decrease of Pw with decreasing 
and increasing effect amplitude w, few trends in the 
values can be easily generalized. From the tabulated data 
we can at least see that for a given combination of q, p,, 
and w the value of P^ is lowest for the d — 1 single large 
caustic case. Such a caustic also has the lowest number of 
protruding cusps when averaged over projection angles. 

* For any orientation except for angles perfectly aligned with 
lines connecting the caustic components. 



To get a better understanding of the dependence on 
lens geometry, we include in Figure [9] plots of the angle- 
averaged widths (Iw), (Icc), and the relative probability 
Pw as a function of lens-component separation d for the 
three mass ratios q and the larger source (p, = 0.02). 
The angle-averaged caustic-crossing-region width (Icc) is 
highest in the single-caustic regime, where it approaches 
or exceeds 1 Einstein radius, and drops to the limits men- 
tioned above after the caustic splits into s everal com- 
pone nts (see also iMao fc Paczvhskil Il991t iNight et"aLl 
|2007( ). The contour widths (1^,) for decreasing w are 
progressively higher and peak later after the splitting, 
but generally follow a similar trend. For small d these 
widths converge to the corresponding contour widths of a 
single-point-mass lens, because the contours surrounding 
the small three-cusp caustics shrink to a point together 
with the caustics. For large d the widths converge to 
double the respective value. 

The right panels in Figure [5] illustrate the change of 
the relative probability P^ with d, thus clarifying the in- 
terpretation of Table [H We first note that the character 
of the curves is the same for all three values of q. In the 
single-caustic regime the values of P^ are low. Still, for 
example Po.oi does not drop below 1.1 for any of the mass 
ratios. P^ increases rapidly to a peak when d moves out 
of the single-caustic interval, and progresses to asymp- 
totic values in the d ^ 1 and d ^ 1 regimes. In the 
wide-binary limit the value of P^ is simply given by the 
corresponding single-point-mass contour width expressed 
in units of source diameter. For the three included con- 
tours and K = limb darkening we get Pq.i = 1-2, 
-Po.oi = 3.4, and Po.ooi ~ H- In the close-binary limit 
the relative probabilities are one third of the correspond- 
ing wide-Hmit values, i.e., Pq.i = 0.42, Pq.oi — 1-2, and 
Po.ooi = 3.6. The quoted values are valid for all three 
mass ratios and both source sizes^. In fact, they re- 
main constant for any source sufficiently smaller than the 
Einstein radius of the lower-mass component of the bi- 
nary. More specifically, from single-point-mass extended- 
source computations we find that for the 1% effect this 
requires p* < 0.2^g/(g -f 1), while for the 0.1% and 10% 
effects the factor before the square root is approximately 
0.05 and 0.5, respectively. 

The character of the {lto){d) and Pw{d) curves for the 

^ The curves for w = 0.001 can be expected to turn over for 
(i < 0.1, the lowest computed separation. 
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Fig. 9. — Angle-averaged projected widths (l^) of regions with absolute extended-source effect |(5cx| > w (left column) and relative 
probabilities P^, (right column) of achieving an |<5ox| > f effect as a function of lens-component separation d for three mass ratios (as 
marlted in each panel) and a p, = 0.02 source. Curves for w = 0.001, 0.01, and 0.1 are marlsed by dotted, short-dashed, and long-dashed 
lines, respectively. Solid lines in the left column mark the angle-averaged projected width of the caustic-crossing region (/cc). The relative 
probabilities are computed in comparison with the probability of caustic crossing, = {lw)/{lcc)- The values of d used in the 
computations are marked by crosses along the w = 0.01 curves. 



small source (p* = 0.002) is similar to the large-source 
results, except for reaching generally lower values. It 
is worth noting here that the scaling of equation (|14p 
with implies that for instance (lo.oi) for a source with 
p* = 0.02 should be the same as (/o.oooi) for a source 
with p* = 0.002, keeping q and d fixed. This can also be 
noticed in Figures [4] and [5l The scaling doesn't hold for 
(lo.i), because the corresponding contour lies closer to 
the caustic than the domain of validity of equation (fn]) . 
In either case Pw will include an additional variation due 
to the increase of (Zee) with p*. 

In order to condense the information from Figure [51 
we average the value of the relative probability over the 
main range of separations in observed two-point-mass mi- 
crolensing events, d G [0.5,3.0]. The obtained values 
of for both source sizes are listed in the last three 
columns of Table [TJ For the larger source we see that en- 
countering a > 1% extended-source effect is on average 
1.43-1.63 times more probable than crossing the caustic, 
the values increasing from g=ltoq = l/9. For the 
smaller source we get a factor of 1.14-1.24 increase in 
probability. It is clear from Figure [5] that these average 
values are particularly sensitive to setting the upper end- 
point of the averaging interval in d, at which Pq.oi keeps 
on increasing. 

5. CHROMATICITY 

Equation (fT4| demonstrated that the amplification de- 
pends on the limb-darkening profile even away from the 
caustic, albeit fairly weakly. In general, the variation 
of stellar limb-darkening profiles with wavelength intro- 
duces a chromatic effect in microlensing. In this section 
we explore the dependence of the microlensing amplifi- 
cation on the limb-darkening profile using the two-term 



PCA model given by equation ([9]). 

By introducing a model of this form in equation ((7]) we 
obtain a limb-darkening-parameter n dependence of the 
extended-source amplification A^{yc, k). The rate of its 
change with k is given by the derivative 



9A,(yc, k) _ Jo fi{r)rdr f2{r)rdr 



DAiVcP*) : 



(25) 

where we introduced the geometry-dependent function 
/s„^o(yc+P*r)/2(r)d2r 



^A(yc,p») 



2nJ^f2ir)rdr 
/s„^o(yc + P*r)/i(r)d2 



(26) 
271" /o fiir)rdr 

The integrals in the numerators of the expression are 
computed over the unit disk of the source E5' . The 
function Da (yc , P* ) can be readily interpreted as the am- 
plification difference between hypothetical sources with 
limb darkening given by f2ir) and fi{r), respectively. 
Of course, the function /2 alone does not describe the 
limb-darkening profile of any realistic source. 

More interestingly, for a given two-term limb- 
darkening basis the dependence on individual parame- 
ters is separated, with the factor preceding Z?/i(yc,p*) 
in equation (|25p depending purely on the limb-darkening 
parameter k, and Da{Yc,P*) depending on the source 
position yc and source radius p*. Equation (PS)) implies 
that at any source position yc the amplification changes 
monotonously with the limb-darkening parameter k. The 
sign of the derivative changes along a curve yc = yac(''') 
implicitly defined by setting DA{yac,P*) — 0. The 
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parameter separation obviously persists in higher-order 
derivatives, and is retained even when computing the 
amphfication difference between sources with different 
hmb darkening (or a source observed at different wave- 
lengths), 

A*{yc,Kl) - A,(yc,K2) = h{Kl,K2)DA{yc,P*) , (27) 



with the K-dependent factor 

h{Kl, K2) — (ki — K2) X 

lo hir)rdr f2ir)rdr 

/o [/iM + n,f2ir)]rdrj^[hir) + ^2f2ir)]rdr 



(28) 



We conclude that for a given source size and a given two- 
term limb-darkening model of the same form as the law 
in equation ([9]) there exists a unique achromatic curve 
yac(T) obtained by setting ^'^(yac/o*) = 0. The ampli- 
fication of a source located anywhere along this curve is 
independent of the specific value of its limb-darkening pa- 
rameter. Note that the curve does depend on the source 
size and on the specific functions /i(r), /2(r) of the limb- 
darkening model. For example, it is somewhat different 
for linear limb darkening and for the PC A model. The p» 
dependence disappears away from the caustic, as shown 
below. 

For the PCA model we use, the variation of /i(ki, K2) 
with Ki^2 is dominated by the (ki — K2) proportionality of 
the numerator of equation (|28p . while the K-dependence 
of the denominator has only a weak influence. Hence, 
the largest value is achieved for the largest difference 
in limb-darkening parameters. As shown in Figure 
the extreme physically permitted k values for our model 
are /^p]< — —0.1620 for the most centrally peaked limb- 
darkening profile^", and Kfi ~ 0.0902 for the flattest pro- 
file". 

Based on the preceding we use these extreme values of 
the limb-darkening parameter k to define the microlens- 
ing chromaticity 

S.Uy.) = , (29) 

where we divided equation (|27p by the amplification of 
a source with k = used in § 21 Thus defined, the mi- 
crolensing chromaticity determines the maximum rela- 
tive variation of the microlensing amplification with the 
parameter of a given limb-darkening model. In other 
words, 5chr(yc) is the minimum photometric precision 
necessary for distinguishing between the steepest and 
the flattest limb-darkening proflles of a source located at 
yc. Note that 6chr depends on non-divergent extended- 
source amplifications, and its behavior close to the caus- 
tic is therefore quite different from (5ex- In particular, 
the change of Schr across the caustic is continuous, and 
its value varies even along the caustic. 

In order to understand the variation of 5chi away from 
the caustic analytically, we expand equation (j29|) for a 

Closely resembling the B-band profile of a range of cool K 
and M giants, e.g., the solar-metallicity ATLAS model with T^ff = 
3500 ii", log g = 1.5, vt = 2kms-^. 

Resembling the /-band limb darkening of a range of O dwarfs, 
e.g., the solar-metallicity ATLAS model with T^ff = 48000 A', 
log g = 5, Vt = 2km s~^. 



small source in regions where Aq is analytic, similarly as 
in S 14.21 We obtain, to first order in pi, 



ScAYc) , (30) 



'5chr(yc)~ /l(Kpk,Kfl) X 

f2irydr fi{r)rdr 



where we utilized the approximate expression for the 
extended-source effect ^cx from equation for a source 
with a brightness distribution I{r) — fi{r). The en- 
tire expression preceding Sex in equation (|30p depends 
only on the limb darkening, hence the source-position 
dependence and source-size scaling is exactly the same 
as for Sex- At a distance from the caustic (and its cusps 
in particular), the geometry of (5chr(yc) contours is thus 
exactly the same as that of the extended-source effect 
(5cx(yc) contours. This also means that the achromatic 
curve yac(''") coincides in this region with the 6cx = 
contour, as implied earlier in § 14.21 and thus it also be- 
comes independent of p,. Using our PCA basis and its 
Kpk and Kfl parameters, we obtain —0.197 as the value 
of the factor preceding Sex in equation pO)) . The con- 
tours of 6cx away from the caustic thus coincide with 
contours of a five times lower value of chromaticity, 

<5chr(yc) ~ -O.2(5ox(yc)- 

It is worth noting here that when defining Schr in equa- 
tion ((29ll . one could just as well divide the amplification 
difference for example by the amplification correspond- 
ing to Kpk, Kfl, or the average of these two amplifica- 
tions. In fact, away from the caustic any of these choices 
would lead exactly to the same approximation given by 
equation ((30)) . The advantage of the choice we made is 
apparent after rewriting equation (|29p in the form 



4hr(yc)= /l(Kpk,Kfl) 



7s^, -4o(yc + p.r)/2(r)d^r 
k,,Mr)d'^ 



/s„^o(yc + p,r)/i(r)d2 



(31) 



Such a separation of the limb-darkening-parameter de- 
pendence from the position dependence cannot be 
achieved by the other possible choices of the denomi- 
nator of equation ([29]) . Only with the present choice 
does the full geometry of the contours remain unchanged 
even when comparing the amplification of sources with 
two arbitrary values of n. This reason aside, from the 
preceding analysis we see that the computational differ- 
ences obtained when using other choices would be limited 
to minor variations of effect amplitude in the vicinity of 
the caustic rather than causing substantial changes. 

The relative ampliflcation difference 612 for two arbi- 
trary profiles can be trivially computed from Schr by a 
simple rescaling. 



1512 (yc)^ 



^*(yc,Ki) -^*(yc,K2) 



^*(yc,o) 

_ h{Kl,K2) 
^(Kpk, Kfi) 



'5chr(yc 



(32) 



in view of equations (P7)) and With values ranging 

from -1 to 1, the /i-factor ratio can be used to recompute 
contour values in the subsequent chromaticity plots, as 
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Fig. 10. — Contour plot of the /i,(ki, K2)/'i(Kpki "fl) factor for 
multiplying (5chr according to equation l|32| l in order to get the 
relative amplification difference S12 of two arbitrary PCA limb- 
darkening profiles with parameters ki and K2. More specifically, 
all chromaticity values in the following three figures can be rescaled 
as required in this way. The factor ranges from 1 in the upper 
left corne r to -1 in the lower right, and is numerically given by 
equation II33I I. 

deemed necessary. More specifically, for our PCA basis 
the ratio is given by 

h{K.i,K2) _ - K2) 

/i(Kpk,Kfi) ~ [0. 2676 + 0.1876(Ki-hK2) + 0. 1316 K1K2] 

(33) 

and is plotted in Figure [TOl 

Before proceeding to numerical results, we emphasize 
that all formulae and equations in this section except 
equation p3p are valid not just for the two-point-mass 
lens, but more generally for any lens and any small 
source with a circularly symmetric brightness distribu- 
tion of the form given by equation ([9]). In the case of 
quasar microlensing, however, one should keep in mind 
that much of the chromaticity is caused by the differ- 
ent extent of the emitting region at different wavelengths 
(jWambsganss fc Paczvhskilll991f ). 

When computing (Jchr(yc) or, more generally, when 
computing A,(yc,K) for different values of n, it is suf- 
ficient to separately evaluate fluxes corresponding to /i 
and /2 , and combine them in the numerator and denom- 
inator of equation ([7]) to get the amplification. In this 
way we obtained for Kpk, Hfi, and k = and plot- 
ted (5chr(yc) in Figure [TT] for the nine lens geometries of 
Figure[T]and a source with = 0.02. 

As discussed above, the geometry of the ^chr contours 
is similar to the geometry of S^x ~ ~5(5chr contours and 
the dot-dashed achromatic curve Schr = coincides with 
the zero-effect curve of Figure |31 except in the immediate 
vicinity of the caustic. Staying away from the caustic, 
regions with positive (5chr, where the amplification of a 
source with a peaked limb-darkening profile exceeds that 
of a flat-profiled source, thus coincide with 5cx < re- 
gions of Figure [3l where the point-source amplification is 
higher than the extended-source amplification. Similarly, 
regions with negative Jchr coincide with dcx > regions 
of Figure [3l The explanation of the positive / negative 



region pattern is analogous to the explanation of the Sex 
pattern in § 14. H taking into account that the peaked pro- 
file gives relatively more weight to the source center and 
less to the limb than the flat profile. Just as in Figure [3l 
extended areas sensitive in this case to limb-darkening 
differences can be found connecting the components of 
compound caustics, even though with lower amplitude 
|(5chr|- Not surprisingly, the dependence of these areas on 
d and q is the same as for the Sex areas discussed in § 14.11 

The detailed views presented in Figure [T^] correspond 
to the same two lens geometries as in Figure [5] with 
source sizes = 0.02 (top left and middle panels) and 
p* = 0.002 (top right and bottom panels). The contour 
geometry close to the caustic is clearly different from 
the extended-source effect Sex contours, as seen best in 
the middle panel. Unlike the Sex = curve, the outer 
branch of the achromatic curve does not coincide pre- 
cisely with the caustic and the Schr > band along the 
caustic is substantially broader than its Sex < counter- 
part in Figure [5l especially close to cusps. The most 
striking feature in Figure [T^] is the negative break at 
yd = 0.146 across the positive chromaticity bands in the 
middle panel. A source on a vertical trajectory extended 
through the break would cross the caustic without chang- 
ing the negative sign of its chromaticity. At the break the 
center of the source is highly amplified due to its position 
on the fold caustic. However, the limb of the source is 
amplified even more because a part of it lies along the 
right vertical fold caustic and another part crosses the 
caustic close to its cusp. The less striking dents in the 
positive band along the caustic on both sides of other 
cusps have a similar explanation. 

Regions with maximum positive chromaticity are lo- 
cated just inside of cusps, where the source center 
achieves maximum relative amplification. Regions with 
maximum negative chromaticity are located about one 
source radius outside folds and outside cusps to either 
side of the cusp axis. At these positions the source limb 
achieves maximum relative amplification, as it passes 
through the cusp or its vicinity and has a large part 
aligned with the caustic. 

The insets in the top panels include cuts through the 
plots along the same three source trajectories as in Fig- 
ure[5l selected to demonstrate the behavior of chromatic- 
ity when crossing the caustic at a fold and a cusp, as 
well as its scaling with source size. The length of the 
source paths in the top right panel with p^, = 0.002 is 
adapted to extend the same number of source radii on 
either side of the caustic. For a source approaching the 
caustic along the fold-crossing trajectory A, the negative 
chromaticity peaks as the limb touches the caustic, then 
reverts to a positive maximum before the center exits the 
caustic, and drops back to a major negative peak before 
fully exiting the caustic, after which the chromaticity 
is practically zero. The first two peaks increase slowly 
with decreasing source size, but the third peak doubles 
from Sehr = -0.14 for = 0.02 to Sehi- = -0.28 for 
p* — 0.002. These high values are caused by the steep 
increase of the extended-source amplification when the 
limb just enters the caustic (in the direction opposite to 
A), with the slope given by the relative brightness of the 
limb. The region with \Sehr\ ^ 0.01 is limited to a band 
from 1 outside the caustic to 1.3 p* and 1.4 p* inside 
for the larger and smaller source, respectively. 
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Fig. 11. — Chromaticity (5chr(yc) for a source with radius p* = 0.02 and PCA limb darkening. Panels correspond to the same lens 
configurations as in Figures [T] and [3] Solid contours are plotted for 5(.hr = ±0.001, ±0.01, and ±0.1 (the latter localized close to the caustic 
and hardly visible here); the dot-dashed contour corresponds to the achromatic curve 5chr{yac) = 0. As marked in the color bar, positive 
values are mapped in shades of orange with black contours, negative in shades of blue with white contours. Areas with |<5chr| < 10~* are 
left white. Circles in panel corners illustrate the source size and the scale difference between columns. 



For sources moving along the cusp-exiting trajectories 
B, the first two extremes are higher than in the cor- 
responding A trajectories, reaching Schr = —0.065 and 
Schr = 0.11 for both source sizes. The first negative peak 
occurs 2.6 from the cusp for the larger source, after 
the limb of the source crosses the caustic above and be- 
low the axis. As in the case of the finite-source effect, 
due to the curvature of the caustic for the smaller source 
the peak occurs already 5.3/9* from the cusp. For sim- 
ilar reasons, the positive peak occurs 0.75 inside the 
cusp for the larger source, and 1.6 from the cusp for 
the smaller source. Outside the caustic the curves dif- 
fer from the fold-caustic case A by slowly dropping to- 
wards zero, with the larger source passing Schi = 0.01 at 
a distance of 2.6 p, from the cusp. The smaller source 



passes ^chr = 0.01 at 6 p* from the cusp. The cusp- 
exiting curves subsequently drop into negative values and 
approach asymptotic from below, as indicated by the 
achromatic curve. Inside the caustic the larger source 
crosses Schr — —0.01 at 3.5 p,, and the smaller source at 
7.3 p* from the cusp. 

The symmetric curves of the cusp-grazing trajectories 
C have a positive peak value Jchr — 0.09 with the source 
center directly at the cusp, skirted by negative peaks with 
^chr = —0.16 at positions ±0.95 p* with the limb at the 
cusp, decreasing in amplitude to —0.01 at ±1.7p* from 
the cusp. As in the case of the extended-source effect 
curves, the decay rate is intermediate between cases A 
and B. The chromaticity curves for the two source sizes 
are practically identical. 
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Fig. 12. — Details of chromaticity i5chr(yc) for a d = 1, q = 1 
lens (top panels), and a, d = 2.05, q = 1/9 lens (middle and bottom 
panels) for source radii p, = 0.02 (top left and middle panels) and 
p* = 0.002 (top right an d bottom panels). Meaning of symbols 
and lines as in Figure [TT] with four additional dotted contours at 
±10-2-^ and ±10-1'^. Top left inset: cuts for source tracks marked 
by arrows; top right inset: cuts for the same number of source radii 
along the same tracks. Caustic crossing occurs at time t = 0; t» is 
the source-radius crossing time. Gray horizontal lines in the insets 
mark contour values <5(.hr = ±10~^, ±10""^'^, and ±10"^. 

Figure [13] illustrates the source-size dependence of 
chromaticity in more detail for the three cuts from Fig- 
ure [12] and a single-point- mass lens for comparison. Fol- 
lowing the layout of Figure [HI the graphs for cuts A, C, 
and for the single lens are plotted as a function of time in 
units of the source-radius crossing time t*, while for cut 
B time is shown in units of tB(p*) = (yO*/0.02)~^/'^ t,. 
In each panel the curves with increasing peak height 
correspond to decreasing source size. For cut A, the 
values of = {0.1,0.05,0.02,0.01,0.005,0.002,0.001}; 
for cut B, = {0.1,0.05,0.02,0.01,0.005}; for cut C, 
= {0.1,0.05,0.02,0.01}; for the single lens, p* = 
{2, 1, 0.5, 0.2}. All the S^hr curves converge for p* ^ to 
limit curves, so that curves for smaller source radii would 
be indistinguishable in the plots. Only the fold-exiting 
t/t^ — > 1 part of our cut A curves would still differ — 
here the convergence is the slowest. The limiting value 
of this negative peak can be obtained from equation (I3ip 
if we replace the ratio of the amplified PGA terms by 
the ratio of limb values /2(l)//i(l). The limiting peak 



value for our PCA basis is —1.46, for linear limb darken- 
ing — 1. Just as in the case of the extended-source effect, 
the fastest convergence to a generic chromaticity curve 
is found in the case of the single-point -mass lens (see 
iHevrovskv et all 120001 : lHevrovskvll2008f ). and the limit 
curves in binary- lens cuts B, C, and most of A are valid 
for sources two orders of magnitude smaller than in the 
single-lens case. Reaching the limit during the exit of 
the trailing source limb in cut A would require reducing 
the source size by several more orders of magnitude. As 
discussed previously for the extended-source effect, due 
to the universality of the local lensing behavior the limit 
curves A, B, C are valid for arbitrary fold, parallel cusp, 
and perpendicular cusp crossings, respectively. Only the 
rate of convergence depends on the geometry of the spe- 
cific crossing. 

Returning to Figure [T^] the extent of the region with 
significant chromaticity along the lens axis in the middle 
and bottom panels is even larger than in the top panels. 
For the larger source in the middle panel the 5chr = 0.01 
contour extends 6.4 p* out from the leftmost cusp, but 
only 0.95 p* out from the rightmost cusp. For the smaller 
source the contour extends 13.8 p* to the left and 2.6 p* 
to the right. The asymmetry has the same cause as the 
(5cx contour asymmetry in Figure [5] — the large difference 
in the narrowness parameters K of the cusps, as shown 

in §11X2] 

For an interesting comparison we refer to iHan &: ParkI 
(|2001h . who present similar plots in their Figure 3 of 
the B — I color-change map and its cuts for the same 
caustic detail and a larger source with — 0.1. Their 
blue peaks correspond to our Schr > 0, their red peaks to 
Schr < 0. As expected from the results above, the peaks 
for their cusp-exiting cut are shifted closer to the cusp. 
A color-change map for the whole caustic of this q = I, 
d — 1 lens and the same large source in their Figure 2 can 
be compared with our chromaticity map for a five times 
smaller source in the top central panel of Figure 1111 

To summarize the results of this section, observation- 
ally significant chromaticity occurs within regions sensi- 
tive to the extended nature of the source, but as a gener- 
ally weaker effect it is concentrated closer to the caustic. 
As in the case of the extended-source effect, for smaller 
sources chromaticity is relevant in regions smaller in 
terms of Einstein radii, which nevertheless extend further 
out from cusps in terms of source radii. Positive regions 
reach highest 6chr values just inside cusps, from where 
they extend along the caustic and along the outer cusp 
axes. Negative regions reach a local extreme inwards 
of the positive band along the caustic, but achieve the 
overall highest |(5chr| in a band skirting the caustic from 
outside. These outer high-negative-chromaticity regions 
are also particularly sensitive to the source size, with a 
smaller source permitting higher differential amplifica- 
tion when its limb enters the caustic. The significance 
of these areas is somewhat reduced by the low amplifica- 
tion at the caustic-crossing onset or end of exit. However, 
in current microlensing monitoring projects caustic exits 
are generally observed with better photometry and sam- 
pling than caustic entries, thus giving ample opportunity 
for measuring the limb darkening of the source. 

The chromatic sensitivity of microlensing is even more 
pronounced if we examine the wavelength dependence 
in more detail, that is, if we proceed from photometry 
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Fig. 13. — Source-size dependence of chromaticity (5chr- Panels correspond to the source tracks marked in Figure [121 as noted in 
the top left corner, and to a source-center-crossing single-point-mass lens (bottom right) for comparison. Curves with increasing peak 
height correspond to source radii p« = {0.1,0.05,0.02,0.01,0.005,0.002,0.001} for panel A, p, = {0.1,0.05,0.02,0.01,0.005} for panel B, 
p, = {0.1,0.05,0.02,0.01} for panel C, and p* = {2,1,0.5,0.2} for the single lens. The time intervals, time units, and vertically marked 



times are the same as in Figure [6] Horizontal dotted lines mark 5chr 

to spectroscopy. Although microlensing spectroscopic 
effects are closely related to broad-band chromaticity 
(jHevrovskv et al.1 12000). their exploration in full two- 
point-mass lensing is beyond the scope of this paper. 
Instead of a simple two-term limb-darkening law such 
a study requires the usage of a full model atmosphere 
with its wavelength-dependent center-to-limb specific- 
intensity variation at the necessary spectral resolution. 

6. VALIDITY OF LINEAR-FOLD APPROXIMATION 

The analysis of stati c binary fold-caustic-crossing mi- 
crolen sing events (e.g., lAfonso et al.ll2000l : lAlbrow et al.l 
l2001a^ with extended-source effect s often follows a proce - 
dure similar to the one outlined bv lAlbrow et al.l (|1999bD . 
An analogous approach has been used for more compli- 
cated events such as OGLE-2002-BLG-069 ( Kubas et al.l 
|2005() , which involves the parallax effect, and in two 
early papers on ER OS-BLG-2000-5 (,Albrow et al.ll2001bl : 
lAfonso et al.ll200il ). a cusp-approaching event that man- 
ifests both binary rotation and parallax effects. 

When a source exits the inner region of the binary caus- 
tic, two of the five images disappear. As a point source 
approaches a fold section of the caustic from inside, 
the amplification of the two vanishing images diverges 



±0.01, ±0.1. 

in a square-root singular ity (|Blandford fc Naravanlll986l : 
[Schneider fc Weij3ll986f) . Assuming that the ampHfica- 
tion of the three remaining images changes only linearly 
near the caustic and approximating the caustic by a tan- 
gent at the point of crossing, the approximate total am- 
plification At^^''{t) of an extended source as a function 
of time t is given by 



At 



+ Acc+Uj{t-tcc), (34) 



where is a scaling factor depending on the source size 
and the characteristic rise length of the caustic, and 
CO describe the linear variation of the amplification of 
the three non- vanishing images, and ice is the time when 
the source center crosses the caustic^^. The timescale of 
the caustic crossing is measured by the time from first 
contact to tec, 

At = i*csc0, (35) 

where — p<,t^ is the source-radius crossing time, t-^ 
is the Einstein-radius crossing time, and 4> is the angle 
between the trajectory and the caustic at the point of 

For a fold-caustic entrance it is s ufficient to change the sign 
of the argument of Gj in equation 1134 l l. 
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Fig. 14. — Caustic-crossing details of events with main parameters taken from OGLE-2002-BLG-069 (left column) and EROS-BLG- 
2000-5 (right column). Top panels: caustic-crossing geometry with overall event geometry in insets. Solid line: caustic; dotted line: source 
trajectory with arrow marking the section used in the bottom panels; dashed line: linear fold tangent to the caustic at crossing point. 
Bottom panels: light-curve detail with fits and their relative residuals <5appr- Solid line: exact light curve A»; dotted and dashed lines: 
best-fit fold-approximation light curves A'^^^'^ for fitting methods I and II, respectively. Bottom right inset: wider part of the light curve 
including the following cusp-axis approach. 

crossing. Gj (rj) is a characteristic profile function which 
depends only on the limb darkening of the source but not 
on its size, 



inax(77, — 1) 



\/r\ - 1] 



dr2dri 



(36) 

where H{1 — rj) is the Heaviside step function, and the 
limb-darkening function /(r) is usually normalized to 
unit total flux. For several classes of analytic limb- 
darkening functions equation l[36|) can be expressed in 
terms of elliptical integrals; in more general cases Gj {r]) 
can be computed numerically. lAlbrow et al.l (|1999br ) ex- 
plain the linear- fold approximation and its parameters in 
more detail. 

A fold-caustic crossing with a fixed limb-darkening pro- 
file can thus be fitted by equation ((34)) with five free pa- 



rameters: U* , tec, A t, Acc, and uj. The procedure of 
lAlbrow et al.l (|1999b[ ) then continues with a fit to the 
non-caustic-crossing part of the light curve constrained 
by the pre-fitted parameters to determine the lens ge- 
ometry. Finally, the full data set is used to refine the 
solutions. However, equation (p4|) is still used for the 
caustic-crossing part of the light curve. 

The outlined methodology has several shortcomings 
that can hinder correct reconstruction of the lens geom- 
etry and the retrieval of source parameters. First, as we 
showed in 21 extended-source effects are not limited to 
the immediate vicinity of the caustic curves. Instead, 
affected areas may be substantially larger (see Figure [3] 
and Table [T]). Second, not every assumption of equation 
(j34p is justified under all circumstances. For example, 
the curvature of the caustic becomes important for tra- 
jectories with small (p, and so does the change of Ao 
parallel to the caustic. As shown in § 14.11 and Figure [SI 
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even in perpendicular crossings the source has to be small 
enough for the light-curve shape to be independent of p, . 
Furthermore, the magnification of the three images un- 
affected by caustic crossing cannot be always regarded 
as a linear function of time, for example when the caus- 
tic crossing occurs in the vicinity of a cusp. Close to 
cusps even the inverse-square-root approximation of the 
divergence breaks down. There have been several cases of 
fold-crossing events, in which equation ([M)) turned out to 
be inadequate for the anal ysis. For example , in the case 
of OGLE-1999-BUL-23 Al brow et all (|2001aD found that 
the difference between the approximation given by equa- 
tion (IM]) and ampli f icatio n computed by the method of 
iGould fc Gauchereil (|1997t ) amounted to as much as 4%. 
To cope with this problem they used pre-calculated tabu- 
lated corrections to improve the approximate expression. 

In the following we compute n— point light curves of 
fold-caustic crossings spanning several At around the 
time of caustic crossing. For simplicity, we assume that 
At and tec are known from the event geometry. Thus, we 
are left with a linear fit of the three parameters C/,, Acc-, 
and Lo. There are two natural ways how to perform the 
fit. We may use standard linear least-squares fitting for 
all n points (hereafter method I). Alternatively, realizing 
that Gi{t]) = for 7y > 1, we let such points determine 
Acc and and points with 77 < 1 are then used to fit [/* 
(hereafter method II). 

We illustrate the drawbacks of the linear-fold approx- 
imation on two caustic-crossing events chosen to resem- 
ble actually observed events, OGLE-2002-BLG-069 and 
EROS-BLG-2000-5. The former event was sele c ted as a 
case favorable for the me thod of I Albrow et al.l (|1999b( ). 
based on the analysis by lKubas et al.l ( '20051 ). We used 
their close-binary parameters (g = 0.58, d = 0.46) but 
we omitted annual parallax because it has a negligible ef- 
fect during the short interval of caustic crossing, and our 
objective is to study the method rather than the event 
itself. The source star with radius p* = 0.0048 followed a 
linear trajectory parameterized by its closest distance to 
the center of mass uq = Mm(yc) = 0.016 and the angle 
a = 134.4° between the yi axis and the trajectory (see 
the top left panel of FigurefMlfor a visualization of the ge- 
ometry). We use here linear limb darkening (eq. [8]) and 
set V — 0.6 2 to reproduce the /-band limb-darkening co- 
efficient of iKubas et al.l (|2005f ). The second intersection 
of the trajectory and the caustic occurs under an angle 
(f) ~ 86°, thus making it a nearly perpendicular crossing. 
We plot the relative residuals (5appr of the approximation 
given by equation (j34p from the exact amplification A^ 
computed from equation ([7]) as described in § [Sj 



^appr . 



(37) 



in the bottom left panel of Figure [M) Even in this fa- 
vorable case the residuals may exceed 2% with either 
fitting method, a deviation from the linear-fold approx- 
imation that is observationally significant. Such devia- 
tions would be deleterious especially for limb-darkening 
measurement, taking into account that the chromaticity 
mostly remains below the level of 10%. With method I, 
<5appr peaks when the edge of the source exits the caustic 
(i.e., when the two highly magnified images disappear) 
and the amplitude of the residual keeps growing over 2% 
as the source moves more than 2At away from the caus- 



tic. When using method II, (5appr keeps growing over 2% 
when the source lies more than 2At inside the caustic. 

In the case of the latter eveii t EROS-BLG-2000-5 , 
initial works (ICastro et all 120011: lAlbrow et all l2001bl : 
'Afon so et al.l [20011) concentrated on the second caustic 
crossing, during which spectra had been meas ured. The 
event wa s subsequent l y thor oughly analyzed bv lAn etldl 
(|2002D . [Fields et all (|2003[ ) refined the event parame- 
ters further and compared the measured limb darkening 
to stellar atmosphere models. In our analysis we adopt 
the parameters obtained by (Fields et al.l (|2003D . namely 
d = 1.940, q = 0.75, = 0.004767, a = 73.851151°, but 
we do not incorporate binary rotation and annual paral- 
lax. Fro m the distance of the closest approach to the cusp 
given by Field s et al.l (|2003f) we compute uq = 0.7096, 
which yields (j) ~ 16.6°. The event geometry is shown in 
the top right panel of Figure [TH Based on the reported 
y-band square-root limb darkening, we use a linear limb- 
darkening profile with v = 0.792^^. A comparison of ex- 
act and best-fit approximated light curves is given in the 
bottom right panels of Figure [TH Obviously, the match 
is poor, locally exceeding 30% for both fitting methods, 
and the reason is clear. Due to the curvature of the 
caustic and the small value of 4> the crossing takes longer 
and, most notably, the caustic exit occurs significantly 
later than would be expected from the linear-fold approx- 
imation. Furthermore, just after exiting the caustic the 
source passes in the immediate vicinity of a cusp, there- 
fore the amplification of the three non- vanishing images 
cannot be regarded as a linear function of time. More 
generally, the amplification strongly changes in the di- 
rection parallel to the tangent linear fold. 

From the bottom right plot in Figure [TJ] we see that 
(Jappr peaks around \t — ted — One might try to 

improve the fit by slightly adjusting tec to account for 
the obvious shift which can be seen in the light curve 
in the right column of Figure [Ml We tried to fit tec to- 
gether with ^cc, and u for both events but found 
only a marginal improvement, with the relative residu- 
als staying in the percent level for OGLE-2002-BLG-069 
and tens of percent in the case of EROS-BLG-2000-5. 
Obviously, one could improve the fits further by varying 
also At. However, this would amount to modifying the 
physical parameters of the event in order to maintain in- 
ternal consistency, and such an analysis is beyond the 
scope of this paper. While it might slightly improve the 
fit for OGLE-2002-BLG-069, the complex pattern of the 
large residuals in the case of EROS-BLG-2000-5 practi- 
cally eliminates any hopes for a plausible linear-fold ap- 
proximation to the crossing. More importantly, it is clear 
that releasing the fixed parameters in the demonstrated 
examples leads to a bias in their recovered values, by try- 
ing to account for the residuals caused by the inadequate 
approximation. 

Despite the preceding discussion, the linear-fold ap- 
proximation could still be useful for analyzing caustic- 
crossing events similar to OGLE-2002-BLG-069 in com- 
bination with an appropriate limb-darkening model — 
but for substantially smaller sources. However, any usage 
of the approximation should be justified by a test of its 

^■^ Computed by a least-squares fit to the two-parameter square- 
root profile, although details of the limb darkening have little im- 
pact on the point raised in this section. 
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applicability in the particular case, at least by comparing 
the exact light curve for the obtained parameter values 
with its best fit by the fold model from equation ([M)) . 

7. COMMENTS ON OBSERVED EVENTS 

According to iMao fc Paczvhskil ()1991l ) ~ 10% of aU 
microlensing events in the Galactic bulge should display 
signatures of strong binary lensing, i.e., they should un- 
dergo caustic crossing. In many non-caustic-crossing bi- 
nary events the influence of lens binarity leads merely 
to perturb a.tions of the point-source - point-lens (PSPL ) 
hght curve (|Di Stefano k Pernalll997HNight et aLll2007D . 
Nevertheless, as we showed in § |4l even in non-caustic- 
crossing events light curves may be affected by extended- 
source effects. In this section we first comment on several 
such observed events that do or potentially could illus- 
tr ate points raised in p revious sections. 

lUdalski et al.l ()2005[ ) analyzed the extensive observa- 
tions of the event OGLE-2005-BLG-071, and concluded 
that the deviations from the PSPL light curve are caused 
by the source passing close to three cusps of a caustic in- 
duced by a Jovian-mass companion to the primary lens. 
Here we assume the following parameters for this event: 
q = 0.0071, d = 1.294, uq = 0.0236, and a = 274.23°. 
In Figure [15] we illustrate the dependence of amplifica- 
tion extended-source effect Sc^{t), and chromatic- 
ity Schi{t) on the source radius / 9<.. The preci s ion o f 
the photometric measurements of lUdalski et al.l ( 20051 ) 
was better than 1%. We can see from Figure [iS] that 
such a precision permits the detection of the extended- 
source effect for sources with > 0.0003, and (5chr 
reaches 1% for sources with > 0.001. This does not 
necessarily mean that meaningful limb-darkening mea- 
surement would be possible in this event, but rather 
that one should select a more-or-less appropriate limb- 
darkening p rofile for the mode ling. The new analysis of 
the event bv lDong et al.l ()2008f ). which appeared after our 
manuscript submission, confirms our conclusions. Source 
radii > 0.0009 are ruled out at > 3 ct, the best-fit solu- 
tion yields = (5.5 ± 1.2) x 10~^, and varying the limb 
d arkening had no sig ni ficant effect on the results . 

iJaroszvhskil (l2002f) jjaroszvhski et all (|2004 I2006D . 
and iSkowron et aLr()2007[ ) published models of 53 bi- 
nary microlensing events observed by the OGLE project. 
Although the number of observations of many of the 
weak binary events is not enough to constrain extended- 
source effects or chromaticity, with better sampling some 
of them would be suitable. F or example, trajectories 
of events OGLE-2004-BLG -280 (iJaroszvnski et al[ [2006h 
and OGLE-2002-BLG-099 paroszvhski et al.ll2004f) cross 
the regions between the caustic components with en- 
hanced sensitivity to the extended nature of the source. 
The same is true for MACHO 97-BLG-41 (jAlbrow et al.l 
|2000[) . but in this case the situation is more complicated 
due to the rotation of the binary. 

Turning to fold-caustic-crossing events, in § [B] we 
discussed OGLE-2002-BLG-069 and EROS-BLG-2000-5. 
We showed that even in the case of a nearly perpendic- 
ular caustic crossing well separated from cusps the bias 
introduced by using the linear-fold approximation is ob- 
servationally significant. We performed a similar analy- 
sis for several other binary events with goo d data cover- 
age. The MACHO-98-SMC-1 analysis by lAfonso et all 
(|2000D yielded two viable solutions, either a close bi- 
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Fig. 15. — Light-curve detail of a planetary-mass non-caustic- 
crossing event approaching three cusps, with main parameters 
taken from OGLE-2005-BLG-071. Top panel: amplification A*; 
middle panel: extended-source eff'ect (5cx; bottom panel: chromatic- 
ity (5(,hr- Individual curves correspond to different source radii p*. 
In the top panel curves with highest to lowest peaks correspond to 
a point source (barely distinguishable), p« = 0.0005, 0.001, 0.002, 
0.005, and 0.01. The order is reversed and the point source omitted 
in the lower two panels (with the largest source producing largest 
deviations). The time axis is plotted in arbitrary units; gray hor- 
izontal lines in the middle panel mark ±10% and ±1% effects, in 
the bottom panel a ±1% effect. 

nary (potentially rotating) or a rotating wide binary, and 
one non-viable static wide binary solution. Similarly to 
EROS-BLG-2000-5, there is a cusp approach just after 
the second caustic crossing. Our inspection of the sec- 
ond caustic crossing for the static close and (non- viable) 
wide binary models shows that i^appr stays on the percent 
level with two peaks at \t — tcc\ = reaching as high 
as 8%. For OGLE-19 99-BUL-23 we assume the preferred 
wide binary model of'A lbrow et al.l (|2001a[ ). Again, S^ppr 
stays below 2% during most of the second caustic cross- 
ing with peaks at \t — ted = /^t re aching 4% at most. 
For this event, lAlbrow etall (|2001a[ ) also found the 4% 
discrepancy from the true amplification. As mentioned 
earlier, they retained the analytic approximation of equa- 
tion (|34p and improved it using pre-calculated correc- 
tions. 

The crossing of a cusp precludes the applicability of 
simple approximations such as in the fold case. Hence, 
the analysis of all such events required the more arduous 
computation of the exact light curve, usually by variants 
of the ray-shooting method. This approach has the clear 
benefit of avoiding potential biases introduced by such 
approximations . 

8. DISCUSSION 

In our analysis we neglected blending, which is a signif- 
icant effect in observed microlensing events. In general 
we expect the observed relative amplification excess over 
the point-source amplification to be less prominent. In 
high-amplification regions near the caustic the blending 
effect will be suppressed for all but the most extreme 
contributions of blended light. The situation is more 
subtle on the axis connecting the caustic components. 
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The point-source amplification here is often rather low, 
e.g., Aq ~ 5, so that adding an ^ 80% blend to the 
baseline may lower Sex by a factor of two. There are sev- 
eral reasons why not to worry about blending too much. 
Commonly used image-subtraction techniques automati- 
cally remove all significant (non- variable) blending. Even 
without image subtraction, being an additive constant 
blended flux is easy to fit and eliminate from light curves 
with decent photometry, leaving the possibility of de- 
tecting (5ex at least in the pattern of residuals. Based 
on the recovered event geometry it should be obvious 
where to look for the effect. In addition, improvements 
in crowded-field photometry and especially any potential 
future space-based microlensing surveys will reduce the 
frequency and impact of strong blending for obser vations 
towa rds the Galactic bulge (see Figure 3 of Han et al.l 
[20061) . 

The current approach to observations of gravitational 
microlensing anomalies is optimized for the detection of 
lenses with exoplanets. In this setup two survey teams 
(OGLE and MOA) image the Galactic bulge at most 
several times per night and provide candidates for fol- 
lowup observations. These are carried out by two other 
teams (PLANET and /zFUN) equipped with a network of 
dedicated telescopes. Followup monitoring concentrates 
on short-lived deviations from the PSPL light curve, es- 
pecially in high-magnification events , where the chances 
for d iscovering a planet are highest (|Griest fc Safizadehl 
Il998f l. Unfortunately, some of the events in which the 
departure from the PSPL light curve does not match 
a planetary signal are ignored. However, many of these 
non-planetary events could provide valuable data on stel- 
lar surfaces that are hardly accessible by other means. 
Based on our results in § |4] we suggest that attention 
should be paid also to the parts of the light curve that 
correspond to approaches to any cusp axis, and the re- 
gions between facing cusps in particular. An analysis 
of these parts may constrain the limb darkening of the 
source even in non-caustic-crossing events and thus en- 
large the sample of stars with measured limb darkening. 

The present observational setup can be expected to 
change in the near future. Wide-field survey telescopes 
will monitor all events in their large field of view with 
high cadence. Recently, the MOA group started to 
observe a subs et of their fields very frequently dur- 
ing the night (jSumil l2008f) . P lanned next-gener ation 
ground-based survey telescopes (jGould et al.l 12007 ) will 
continue in the trend of dismantling the alert/foUowup 
modus operandi. Proposed space missions such as the 
Microlensing Planet Finder ([Bennett et al.l [20071 ) would 
provide an almost nonstop flow of high-quality data on 
microlensing events in the bulge. As we have shown in 
§ [H for a giant source the area sensitive to its extended 
nature can be up to several times wider than the caustic- 
crossing area. For main-sequence sources this width in- 
crease can reach several tens of percent. The develop- 
ments mentioned above would give excellent prospects 
for a massive survey of limb darkening of bulge stars 
with two-point-mass microlensing. 

As we demonstrated, modeling of binary microlensing 
events requires a very accurate computational method. 
Poor approximations introduce bias not only in the mea- 
sured limb darkening but also in the recovered binary 
parameters. All of the events in which modeling adopted 



the linear-fold approximation given by equation (|34p 
would benefit from a re-analysis incorporating an accu- 
rate light-curve computation method. This applies in 
particular to events used for source limb-darkening mea- 
surement, as the positions with highest chromaticity are 
largely co ncentrated at arou nd the caustic where the 
method of lAlbrow et al.l (|1999b[ ) tends to produce largest 
deviations. 

9. SUMMARY 

We imp lemented the numerically efficient method of 
IVermaaM (2000, ') for computing the general two-point- 
mass microlensing amplification of a source with an ar- 
bitrary surface-brightness distribution. Using a PCA- 
based model of stellar limb darkening, we explored the 
sensitivity of such lenses to the extended nature of the 
source and to differences in its limb darkening. Tradi- 
tionally, all extended-source effects were expected to be 
limited to a narrow band along the caustic. In contrast, 
we discovered regions of strongly enhanced sensitivity 
along the outer axes of cusps, which may even bridge 
the gap between caustic components of close or wide bi- 
nary lenses. Although for smaller sources the extent of 
the sensitive region along the cusp axis shrinks in terms 
of Einstein radii, we show that in terms of source radii 
it grows, scaling with the inverse cube-root of the source 
radius. 

We derive analytical approximations for the extended- 
source effect and chromaticity of a small source not posi- 
tioned directly on the caustic. These expressions can be 
readily applied not only to the two-point-mass lens, but 
more generally to any other gravitational lens system. 
We use exactly computed light curves to check the ap- 
propriateness of the linear-fold approximation, which is 
often used in the analysis of caustic-crossing events. The 
approximation leaves residuals on the level of a few per- 
cent even in near-ideal events such as OGLE-2002-BLG- 
069, which meet most of the conditions of the approxima- 
tion. It fails badly in events such as EROS-BLG-2000-5, 
which violate the conditions excessively. All microlens- 
ing events in which the approximation was employed for 
measuring the limb darkening of the source star there- 
fore deserve a re-analysis, if only as a consistency check 
in the more favorable cases. Any usage of the linear-fold 
approximation in event analysis should be accompanied 
at least by a simple test demonstrating its applicability 
in the particular situation. 

A study of the size of contours of the extended-source 
effect illustrated that for instance the probability of a 
microlensed bulge giant causing at least a 1% effect is on 
average 1.4-1.6 times higher than the probability of its 
caustic crossing, when averaged over a range of typical 
lens-component separations. Overall, our results bode 
optimistically for the prospects of measuring the limb 
darkening of source stars in a substantially higher num- 
ber of two-point-mass microlensing events. 
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y, y[K] 

Fig. 16. — Geometry of a generic cusp caustic given by equation I IA2I I. Left panel: cusps for six values of K, with K > curves drawn 
bold and A" < thin. Cusps with \K\ = 1, 8, and 0.125 are marked by solid, dashed, and dot-dashed lines, respectively. Right panel: 
amplification contours for lensing by a generic cusp (bold curve). Contour values in units of l/Cbl"^ are spaced by a factor of 2^/^; outside 
the cusp levels range from the outermost 2~^/^ to the innermost 2^ , inside the cusp from the rightmost 2^/"^ to the leftmost 2^ . With axes 
marked in units of K , the right panel is valid for an arbitrary cusp. 



APPENDIX 



AMPLIFICATION OF IMAGES FORMED BY A CUSP 

As shown by [Schneider fc Weifil (|1992f ). the lens equation near a generic cusp caustic placed at the origin may be 
written as 

h 



y\\ =cx|| -1 



-xi 

2 ^ 



(Al) 



where (y||,yj_) are components of the source position parallel and perpendicular to the axis of symmetry of the cusp, 
(x||,xj^) are analogously defined coordinates of the image position, and a, 5, c are constants fulfilling the conditions 
6 7^ 0, c 7^ 0, and — 2ac ^ 0. The equation of the cusp caustic describes a semi-cubical parabola 

(A2) 



where 



K 



27c2 
"863' 



[b^ - 2ac) 



(A3) 



is a parameter determining the orientation of the cusp and its narrowness, as shown in the left panel of Figure 1161 
A cusp with K > \s pointed to the left (towards negative yy), a cusp with if < is pointed to the right (towards 
positive y||). A cusp with a large value of \K\ is narrower than a cusp with a small value of \K\. Without a loss of 
generality we choose for the rest of this section the orientation with K > Q. If this is not fulfilled automatically by 
equations (|Aip . it can be readily achieved by rotating the coordinate systems y and x by 180°, which corresponds to 
changing the sign of parameter b. 

Inversion of equations (jAl[) leads to a cubic equation for x j ^, and the ampli f icatio n of each image can be expressed 
as a function of the image position (jSchneider fc Weij3ll992f ). IZakharovl (Il995l . fl999f ) presented an elegant alternative 
method for obtaining image amplifications as a function of source position, without explicitly solving the lens equa- 
tion. From general properties of the roots of the image equation Zakharov derived a simple cubic equation for the 
amplification, the roots of which are the amplifications of the images. We reproduce the main results here because of 
a few typos appearing in the original papers. A point source satisfying y\ < K^^ y^ lies inside the cusp and has three 
images with amplifications 



yii 



cosi 



arcsm , 



yff 



2^(j - 1) 



J = 1,2,3 



(A4) 



Note that these amplifications include the appropriate sign corresponding to the parity of the image. The j ~ 1 image 
parity is equal to the sign of b, while the j = 2 and j — 3 images both have the opposite parity. The total amplification 



Extended Source in Binary Microlensing 



25 



of a source inside the cusp = X]j=i l^'"*^! = 2 hence we have 



yii 



- arcsm^ 



(A5) 



A point source lying outside the cusp satisfies > K y| and has a single image amplified by 



yl - y3 + y^ 



1/3 



yl-K-^yl-yA 



1/3 



2h 



(A6) 



This image has the same parity as the j — 2 and j — 3 images, that is, opposite to the sign of b. The total amplification 
of a source outside the cusp thus is 

^'°*(y||,y^) = |A(°)(y||,y^)|. (A7) 

In the right panel of Figure [16] we plotted contours of constant total amplification (y) with levels given in units 
of \K b 1^^. By marking both axes in units of K the plot is valid for an arbitrary K > cusp. 
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